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ABSTRACT 

The characteristics of computer based instrumentation, or 
smart instruments, are investigated. Computer based 
instruments are distinguished by their ability to include a 
more complex model of the physical processes influencing the 
desired measurement than is possible with conventional 
instrument. Smart instrumentation is described with 
emphasis on illustrating its ability to improve data 
collection, storage, display, and evaluation. The 
application considerations of redundancy, consistency, 
noise, and filtering are also addressed. 

As an application example, a smart instrument for measuring 
steam generator water level in a pressurized water reactor 
is designed. A model, accounting for the important 
processes affecting level measurement is developed. An 
error exists in the computation the fluid shear stresses, 
but the model calculations remain illustrative of those 
pertinent to smart instrument design. The model is stable 
for both steady state and transient conditions, but there 
are restrictions on the rate of the transient. Simulated 
level data is used to compare a simplified level instrument 
with the smart level instrument. The smart instrument is 
more accurate, but not by more than one percent. Methods 
that could be used by a smart instrument to recover from 
operation outside its model assumptions are discussed. 
Finally, recommendations are made for future work. 
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Introduction 



Chapter 1 



The objective of this research is to suggest a method for 
the application of microcomputer based instrumentation to 
nuclear power plants. The specific instrument illustrated 
is for a steam generator level measurement, but the basic 
steps of physical process selection, variable 
identification, model creation, and integration of the 
process model with signal validation are applicable to any 
important plant parameter. A procedure to implement a smart 
instrument is suggested. A model was developed to include 
important physical features of the level instrument process. 
This model was simplified and used as a basis for a computer 
based steam generator level instrument . Although the model 
is based on fundamental properties and relationships, its 
accuracy and the limits of its application cannot be 
established fully until comparisons can be made with 
experimental data. 

1 . 1 Purpose of Simulation 

The optimum method for evaluation of the proposed steam 
generator level instrument would be to compare the outputs 
of a conventional level measurement device and the 
measurements of a microprocessor-based level instrument 
given identical inputs from an actual steam generator. 
However, since the concept is still being developed, the 
more prudent approach is to perform extensive lab tests 
prior to implementation of this smart instrument on an 
actual power plant. To evaluate the potential effectiveness 
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of the proposed smart level instrument, typical steam 
generator level pressure and power relations are established 
to provide simulated "measured" data for the physical input 
variables that the smart level instrument would require and 
for a simulated "known" level signal for the basis of 
comparison. While it is somewhat artificial to base 
measurement results upon simulated input data, this 
procedure does provide a qualitative and quantitative sense 
for the capabilities that a smart instrument would have and 
serves as a first step for future research. 

1 . 2 Smart Instrumentation 

In the context of this thesis, a smart instrument is one 
that includes a detailed model of the physical process being 
measured to enhance both the reliability and accuracy of the 
instrument. In addition, the model can give the instrument 
the capability to infer quantities that are not actually 
measured. For example, a temperature sensor could be used 
to control a heating coil for bringing a solution in a 
beaker to a desired temperature. A smart temperature sensor 
could detect the absence of the solution from the beaker and 
protect the heating coil by noting whether the rate of the 
temperature increase just after energizing the heater was 
consistent with that expected when the solution was present, 
and de-energize the heating coil appropriately. 



1 . 3 Background 

Accurate indication and measurement of plant variables are 
required for safe operation of nuclear reactors . 

Information recorded by instruments is used to advise plant 
supervisors about the operability and reliable operation of 
nearly all plant components, to provide inputs to protection 
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systems, and to furnish input to automatic control systems. 
In the past, instrumentation was not sophisticated enough to 
enhance accuracy by using information about the process 
being modeled; often attempts were made to design 
instrumentation resistant to sources of signal degradation, 
but the design was usually based on worst case or very 
simplified assumptions about sources of measurement error. 

In addition, conservative estimates are generally made to 
account for the tendency of the transducer and associated 
circuitry to exhibit characteristics varying with time so 
that regular calibrations could be done prior to receiving 
erroneous readings. Recent attempts to integrate 
microcomputers into plant instrumentation and to employ 
signal validation techniques have illustrated the potential 
benefits of more sophisticated instruments (references SI 
and M2) . Signal validation increases information 
reliability by combining information from redundant sensors 
and a detailed knowledge of the dynamics of plant systems. 
Detection of faulty sensor operation, a highly desirable 
function, is an important component of signal validation 
that can be accomplished by high speed comparison of 
redundant sensors in conjunction with knowledge of prior 
sensor data. 

The premise of this research is that microcomputers, when 
included as integral components of plant sensors, can 
perform the functions of signal measurement, validation, and 
fault detection to enhance instrument reliability and 
accuracy. Using microcomputers within each sensor enables 
one to program into the instrument an accurate model of the 
physically pertinent variables and their relationships which 
can enhance the accuracy of the instrument as well as 
provide a localized determination of the quality of data 
presented, also known as distributed processing. This 
information can be used to alert the operator to the onset 
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of conditions that could produce erroneous readings as well 
as extending the time between required instrument 
calibrations. In addition, redundancy and consistency 
checking of measured data can be done locally by a group of 
smart instruments . 



1 . 4 Motivation 

Accurate steam generator water level information is 
essential for safe operation of pressurized water reactors 
(PWR) . The steam generators are the source of the steam 
used to drive the power turbines and are the heat sink for 
the primary coolant that flows through the reactor. 
Maintaining steam generator water level below the upper 
limit of the designated operating band is necessary to 
prevent carryover of liquid droplets to the turbines, which 
could cause severe turbine blade erosion and possibly 
tremendous damage to the turbine itself. Steam generator 
water level must not be allowed to drop below the lower 
limit of the designed band since this could uncover the 
generator U-tubes, resulting in inadequate heat removal from 
the reactor and significant thermal stresses to the U-tubes 
themselves, possibly leading to their premature failure. 
Clearly, it is desirable to enhance the reliability of the 
level information provided to plant operators and to the 
automatic control systems responsible for maintaining 
generator water level within normal bounds . Using a 
computer to record plant data is not a new concept, but the 
smart instrument would be able to interpret the measured 
data as well and by placing the computing resource as close 
as possible to the variables being measured, higher level 
computers utilized in the plant can be freed to perform 
integrating roles on the basis of input information of 
improved quality . 
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Chapter 2 



Steam Generator Description 



Each primary loop in a Pressurized Water Reactor contains a 
steam generator, important and very large (typically 20 
meters high, 4 meters in diameter) components which are used 
to transfer heat from a primary coolant, flowing through the 
reactor core at high pressure, to a secondary coolant, which 
does not flow through the core, to produce the steam 
necessary to drive turbo generators. 

2 . 1 Steam Generator Internal Description 

There are three basic steam generator designs: vertical 

U-tube, once-through, and the horizontal U-tube. Since the 
vertical U-tube generator (Figure 1) is the most widely 
used, only its operation and level sensing configuration 
will be discussed here. The steam generator is a heat 
exchanger that enables the primary coolant to boil the 
secondary coolant without allowing the two fluids to mix. 

The hot, primary coolant just exited from the reactor vessel 
enters the steam generator at the inlet plenum, located at 
the lower end of the steam generator, passes upward then 
downward through the U-tubes, enters the exit plenum and 
returns to the cold leg reactor piping. Mixing of the inlet 
and outlet plenums is prevented by an isolation plate 
separating the two . Condensed steam, known as condensate or 
feedwater, enters the steam generator via an annular piping 
arrangement inside the generator called the feedwater inlet 
ring. Cold feedwater is prevented from impinging upon the 
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tubes by a tube bundle wrapper which forces the feed water 
through the downcomer region, so the feed water can be 
pre-heated by the hot water falling from the moisture 
separators. The feedwater flows up through the U-tube 




region where it is heated and transformed to steam because 
of the lower pressure of the secondary loop. The two phase 
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mixture passes upward from the U-tube region through the 
moisture separators, which allow the steam to pass but cause 
virtually all of the entrained liquid to fall back within 
the steam generator and mix with incoming feedwater, as 
mentioned above. 




Figure 2 : Steam Generator Level Sensing System 

2 . 2 Level Sensing 

Figure 2 shows a schematic representation of a typical steam 
generator level remote sensing system. It is analogous to 
the method used to measure core water level in a boiling 
water reactor (BWR) and the method used to measure the 
pressurizer level in a pressurized water reactor (PWR) 
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plant . The system depicted is sometimes called a cold 
reference leg system because the reference leg piping is not 
heated or insulated so it assumes the temperature of the 
ambient air in the reactor building. The downcomer level, 
the level in the area between the tube sheet wrapper and the 
generator vessel, is the variable actually measured and is 
considered the generator water level. That level is more 
easily measured and the two phase boiling process in the 
tube bundle region makes the level there ill-defined. As 
shown in Figure 2, an upper tap in the steam space is 
connected to a condensing pot and a reference leg. A lower 
tap in the liquid space is connected to a variable leg. The 
pressure difference, a P, between the reference and variable 
legs is measured by a differential pressure cell. The 
liquid level L may be calculated under the following 
assumptions : 

(1) the fluid densities in the reference and variable 
legs are known, 

(2) the liquid and vapor are separated, with liquid at 
the bottom and vapor at the top, 

(3) the reference leg is completely filled with liquid 
to a known height, 

(4) the density of the reference leg is constant, and 

(5) there are no velocity dependent contributions to 
the pressure differences. 

Figure 3 shows a schematic of the differential pressure cell 
arrangement, simplified to show just the details associated 
with the variable and reference legs. Based on figure 3, the 
differential pressure is approximately: 
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Figure 3: Differential Pressure Cell 



P A -P B = *P= (Pc-Psg) + P r Hg-p dc Lg-p v g(H-L) (1) 



and so the level is: 



( PcZEsg ) 

( P de — P v ^ 3 



+ 



( Pr-Pv ' 

\ Pde - Pvy 



H 



( 2 ) 



where : p c 4 

P sg * 

Pr A 

Pdc A 

Pv A 

g 4 

H 4 



L 4 



pressure in condensing pot; 
pressure in steam generator; 
density of fluid in reference leg; 
density of fluid in the downcomer; 
density of vapor in the steam generator; 
acceleration due to gravity; 
height of water in the reference leg; 
and 

downcomer level . 



Note that the level calculated by equation (2) is 
approximate since the height difference between the level of 
liquid in the condensing pot and the reference leg tap on 
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the steam generator is neglected. Otherwise, the numerator 
of the second term on the right hand side of equation (2) 
would be (p z H- p v H tt ) where H tt is the vertical distance 

between the reference and variable leg taps on the side of 
the steam generator. Since p r is so much larger than p v , 
this approximation introduces little error. 

2 . 3 Limitations of Conventional Level Sensing 

Generally, the assumptions used for level measurements and 
noted above are valid and equation (1) provides a good 
estimate of steam generator level during normal operations. 
However, severe transients involving rapid depressurizations 
such as a generator steam leak could lower the pressure 
above the reference leg enough to cause flashing in the 
reference leg, violating assumption 3. Variations in the 
temperature of the reactor building will affect the rate of 
condensation in the condensing pot and the average density 
of the liquid in the reference leg. Generally, reactor 
building temperature is moderately related to time at power 
the power level of the reactor, and the outdoor temperature. 
In addition, a casualty such as a steam leak from the 
generator or a primary coolant leak could alter reactor 
building temperature enough to make assumption 4 invalid. 

The diameter of the piping chosen for the condensing pot 
drain line is also an important consideration. If the 
piping is not large enough to accommodate the flow of 
condensing liquid from the condensing pot during worst case 
conditions, the pipe may become filled with liquid, 
resulting in a decrease in effective reference leg height as 
the water level in the condensing pot rises above the drain 
line piping connection to the condensing pot. This 
condition violates assumption 4 . There have been 
well-studied examples of this phenomenon actually occurring. 
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Finally, the steam that flows from the steam generator to 
the condensing pot and condenses does so because the 
condensing pot is uninsulated, but it would not be accurate 
to describe the condensation as taking place at reactor 
building ambient temperature. Contact between the inner 
surface of the condensing pot material and the steam from 
the steam generator no doubt elevates the average 
temperature of the condensing pot surface above reactor 
building ambient temperature. Still, the portions of the 
reference leg piping in the vicinity of the differential 
pressure cell are most nearly at ambient for the reactor 
building. Clearly, assumption 4 is not really accurate and 
the density of the reference leg varies along the length of 
the reference leg piping. 
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Chapter 3 



Smart Instruments 



A smart instrument is an electronic device, possessing 
memory and controlled by a microcomputer, used to measure 
some plant variable. Its computational and decision making 
capabilities, the salient features distinguishing smart from 
conventional instruments, enable the designer to integrate 
within the instrument a model of the physical process 
affecting the parameter being measured. To implement a 
smart instrument, the designer must choose an appropriate 
model and consider the roles of measurement consistency, 
data redundancy, and filtering. Smart instruments have many 
potential advantages over conventional instrumentation, but 
there are some tradeoffs to be considered as well 




Figure 4 : Smart Instrument Block Diagram 



3.1 Model 

Figure 4 shows a block diagram of a typical smart 
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instrument. It takes its inputs in the form of physical 
variables from some actual process and translates these 
natural phenomena into electrical variables like voltage or 
current. Conventional instruments process a transducer 
output signal in accordance with pre-determined assumptions 
selected by the designer. The processed transducer data is 
routed to a meter face and possible also connected to a data 
recording computer. While this has been a generally 
satisfactory method for presenting plant variable 
information to operators, it has two drawbacks. First, the 
burden of evaluating the consistency and accuracy of the 
data presented is entirely upon the operator. This context 
evaluation burden is not onerous when the total number of 
variables that the operator has to monitor are small, but it 
becomes considerably more difficult in the control room of a 
commercial power reactor having thousands of instruments and 
potentially hundreds of alarms to monitor. 

Detecting instrument or transducer failure is especially 
complicated when the instrument continues to display a 
normal reading although its ability to respond at the 
extremes of its range might be severely degraded. The 
second limitation of conventional instruments is that they 
are only accurate as long as the design assumptions remain 
valid. The microcomputer controlling the smart instrument 
enables the designer to program into the instrument an 
"understanding" of the numerical characteristics of the 
system being measured. The complexity and accuracy of the 
model chosen for inclusion into the smart instrument is an 
important consideration. The primary benefit of a good 
model is that it enables the instrument itself to monitor 
measurement results and alert the operator if the data 
received from the sensors is not consistent with either: 
redundant transducers measuring the same variables or with 
other sensor information that can be fed to the smart 
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instrument. Judicious choice of the parameters to be 
monitored by the instrument and included into the 
computational model is essential. Selecting too few 
parameters and an oversimplified model risks realizing only 
a marginal improvement over existing instrumentation while 
an exceptionally complicated model that requires numerous 
interconnections with other sensors would require an 
enormously powerful computation engine that might dwarf the 
size of the original instrument, posing considerable 
debugging and robustness challenges. 



3 . 2 Data Collection 

It is important to establish some confidence that the data 
collected by the smart instrument is reliable. The 
important issues are: consistency of the detected signals 

with previous signals from the same transducer, comparison 
of the detected signals with current redundant signals, and 
filtering the detected signal to reduce the influence of 
electronic or process noise. 

3.2.1 Consistency 

It is desirable that the smart instrument be capable of 
checking each measured data value for consistency with other 
variables. Consistency checking of instrument results is a 
key supervisory task performed by the plant operators and it 
is plausible that much data supervision could be performed 
by computer controlled instruments . The instrument should 
have memory of past measurement results so it can compare 
the current measurement with recent values of the same 
measurement or some average of the recent values . The smart 
instrument can check the current signal trend with the trend 
exhibited by other variables. The smart level instrument 
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could be programmed to regard a deviation of the current 
steam generator level datum from the average of prior 
readings in excess of some threshold as normal if steam 
generator pressure is changing at some rate and abnormal if 
it is not. Incorporating a process model into the level 
instrument gives it the ability to determine if a particular 
measurement result is consistent with other plant variables. 
Since a smart instrument is envisioned as an aid to the 
operator, it should not require additional alarms or require 
significant retraining to use (reference SI) . The results 
of the consistency checks it performs could be presented to 
the operator by means of a small warning message on the 
instrument meter face or cathode ray tube display as well as 
being fed to a plant data logging computer. 

3.2.2 Redundancy 

Redundancy must be designed into any instrument that 
monitors an important reactor plant variable (reference M2) . 
It would be a serious shortcoming for an instrument 
providing display of and control system inputs based on a 
plant variable measured by a single sensor since it could 
lead to undesirable operator or control system response 
either during normal operations (e.g. failure to act on a 
malfunction because it is not indicated by the instrument) 
or in casualty situations (e.g. masking the casualty or 
giving inconsistent output) . 

The smart level instrument should have at least 2 
independent level sensing signals as inputs. It is already 
a common procedure in power plant design to have 4 
independent electric level sensing channels per steam 
generator so this is not an overly restrictive requirement. 
The smart level instrument would record and compare the two 
level signals at each sampling interval. Level data points 



26 



could then be averaged over a suitable number (selected to 
establish a desired confidence interval) of sampling 
intervals. In addition, the smart instrument could be 
programmed with an allowable deviation from the average of 
prior signals, either from prior knowledge or using 
information recorded in real time, so that sensor data 
points that deviated in excess of the allowable amount from 
would not be included in the average. If a certain number 
of consecutive signal samples are rejected or some threshold 
of total signal samples over a certain time interval is 
exceeded, the operator could be alerted so the affected 
instrument channel could be checked. 

3.2.3 Filtering 

The instrument designer must anticipate that the sampled 
data presented to the smart instrument will be subject to 
noise degradation. The noise would result from interference 
from other electrical signals, process noise, and 
quantization noise. 

Radiation shielding considerations for the reactor and 
reactor building dictate that there be a minimum of 
penetrations in the secondary reactor shield to enhance 
shield structural integrity and reduce the possibility of 
radiation streaming. Minimizing the number of shield 
penetrations necessitates routing multiple electrical signal 
cables in bundles through a limited number of shield 
penetrations. If the smart instrument is located outside 
the secondary shield, it is a virtual certainty that cables 
connecting the smart instrument to its transducers will be 
placed in proximity to higher voltage and or frequency 
signal cables for other systems so that there will be 
crosstalk from these other signals that will reach the smart 
instrument input. There will probably also be a significant 
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noise component at 60 Hertz due to power line interference. 



Process noise is composed of random signal variations 
introduced at the transducer input. The random component of 
the signal is produced by minor local variations in the 
measured values and the statistical nature of the physical 
interactions between the measured variables and the 
transducer. In addition, the smart instrument can be 
programmed to evaluate sensor performance by interpreting 
the character of the process noise. The instrument could 
either be pre-programmed with maximum and minimum tolerable 
values for noise variations or could determine these limits 
through a weighted average of prior data. This expected 
noise deviation data could be used to alert the operator to 
an abnormal sensor noise signature and thus detect otherwise 
undetectable sensor failures. 

Quantization noise results from the inability of a digital 
computer to accept data in the continuous, analog form 
produced by the transducer. Before the smart instrument can 
process the sampled signal, it must first be converted to a 
digital representation. The resolution or accuracy of the 
analog to digital conversion is related to the number of 
binary digits, bits, used to represent the digital data. 
While quantization noise can be reduced by using an analog 
to digital (A/D) converter with more bits (at the cost of 
storing and transmitting more data) , it is an important 
consideration for the designer, especially since the digital 
signal values will be used in various computations of model 
parameters . 

To compensate for process noise and interference noise 
degradation of the sampled data, the smart instrument could 
be programmed to filter the incoming data digitally prior to 
processing in the model. It could also accept and record 
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the raw, unfiltered, digital data so it could be compared 
with the filtered data if necessary for conducting 
performance diagnostics. One could select a filter to 
remove a range of frequencies or reject a specific frequency 
where the majority of the noise power might be concentrated. 
A moving average filter would be useful for getting a good 
estimate of the mean signal value in the presence of noise. 

One other benefit of using a digital filter is that the 
filter characteristics could be easily changed without 
requiring a redesign of the instrument or the filter could 
be adaptive, possibly using a matched filter, to adjust the 
filter parameters in real time as the information is being 
received . 

3 . 3 Capabilities and Limitations 
A smart instrument's capabilities include: 

'Enhanced Reliability: a smart instrument can provide 

more reliable information than conventional sensors because 
it can include some extra variables in its model that enable 
it to provide reliable indication over a broader range of 
situations, needs fewer assumptions, and can correct 
transducer drift if the designer has the ability to model 
it . 



•Calibration Reduction: smart instruments provide a 

real time means for objectively assessing the quality and 
accuracy of data presented, thus replacing strictly calendar 
based calibration requirements. 



•Variable Inference: in some situations, a smart 

instrument can infer plant variables that are not directly 
measurable, but are desired for increasing model accuracy. 
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assessing instrument reliability, or for providing extra 
operational information. 

•Self Diagnosis: a smart instrument can be programmed 

with built in test (BIT) routines to perform diagnostics on 
its major components and alert the operator to specific 
failures, thus reducing repair time. 

•Sensor Communication: smart sensors can query other 

sensors for information on a common communication bus 
without requiring dedicated instrument to instrument 
connections. A control circuit could periodically poll 
sensors for data. Sensors could record the time when a 
certain reading or alarm occurs to be used for fault 
isolation. For example, pressurized water power plants 
usually have salinity cells at various locations in the 
condensate stream of the secondary system to detect the 
introduction of impurities into the condensate. 
Unfortunately, a large sea water leak upstream of several 
salinity cells will usually cause all the downstream cells 
to alarm nearly simultaneously so it is difficult to specify 
rapidly the source of the contamination, thus requiring more 
time to find and correct the source of the casualty. If the 
salinity sensors recorded the times when their alarms 
occurred, the operators might be able to identify more 
easily the source of the contamination based on the recorded 
alarm times. In addition, a communications channel between 
instruments could enable a control unit or protection 
circuit to query the instruments that provide it with input 
data, either to perform continuity checks or brief on-line 
diagnostics. Depending on the sampling interval and the 
hardware used, there could be sufficient time between data 
samples to do these checks . 

The limitations of smart instruments include: 
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•Model: the performance benefits discussed above are 

limited by the quality and accuracy of the model used by the 
smart instrument. The model must be chosen carefully so 
that it is robust and not subject to oscillatory or unstable 
behavior . 

•Model Complexity: the designer does not have 

unlimited flexibility when selecting model complexity and 
variables to include. The more complicated the model, the 
faster the central processing unit needed and the more 
memory that will be required. In addition, robustness tends 
to be a more serious issue for more complicated models . 

•Fault Tolerance: the smart instrument must be capable 

of withstanding power failures and a potentially harsh 
operating environment. Its data presentation algorithm and 
model implementation must be capable of recovering from 
unexpectedly large and small inputs without suffering 
failure. The designer might decide to use some kind of 
battery backup for the instrument, especially if it does 
some of its processing based on real time variable history. 

•Interface: computers and computer controlled devices 

generally exhibit the tendency to produce and display data 
that has the appearance of truth even when receiving faulty 
inputs. It is important that the data be presented to the 
user in a form that enables him to make a determination of 
the quality of the data. 

Conclusions 



Because it is controlled by a machine that is capable of 
making complicated decisions based on input data and has the 
ability to adapt its performance to the measured data during 
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operation, the smart instrument is capable of considerably 
more sophisticated data analysis and manipulation than a 
conventional instrument. 
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Chapter 4 



Dynamics of a Steam Generator Level 

Sensor 



The model is developed in this research to relate the 
pressure in the condensing pot to that in the steam 
generator. The model permits calculation of the superficial 
vapor and liquid velocities and the liquid height in the 
drain pipe for the two phase flow from the condensing pot to 
the steam generator. None of these quantities is normally 
measured, but a level instrument that could accurately 
estimate their value dynamically could alert the operator to 
the onset of conditions that could lead to inaccurate level 
information . 

4 . 1 System Description 

The purpose of the condensing pot is to ensure the reference 
leg of the level instrument is always filled to a constant 
level. During normal operation, steam flows from the 
generator steam space to the condensing pot via a piping 
connection, called the drain line or tube in this paper. In 
the condensing pot, the steam condenses because the 
condensing pot is not insulated (as noted in earlier 
descriptions) . Some of the condensed liquid keeps the 
reference leg full, but the remainder spills over from the 
condensing pot and into the drain pipe, returning to the 
steam generator. The condensing pot is modeled as cylinder 
with its axis horizontal having the drain pipe connected to 
one of the circular bases. 
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A design aim for the drain pipe is that it be large enough 
to accommodate the counter-current vapor-liquid flow and 
remain smooth-stratified or wavy-stratif ied under most 
normal conditions. However, some combinations of drain pipe 
diameters and phase flow rates in the drain pipe can lead to 
a flow transition from wavy— stratified to either 
intermittent , in which waves at the surface are picked up to 
form a frothy slug which is propagated along the channel, or 
annular-dispersed, the formation of a vapor core with a 
liquid film around the periphery of the pipe. 

Rearranging equation (1) produces: 

p a~ p b= (Pc-P 3g ) + (Pr - Pv)9 rif -(Pdc-Pv)9 rI ' O) 

As noted earlier, erroneous level indication could result 
from failure of the reference leg height to remain constant. 
As shown in figure 5, if the height of liquid in the 
reference leg rises above the design level (which would most 
likely a result of flooding in the drain pipe) , the 
differential pressure transmitter would sense a decreased 
differential pressure since the effective H is actually 
lower. Since H is assumed by the instrument designer to 
remain constant under all circumstances, the instrument 
would interpret this change as a decrease in level . The 
effect of this condition on the instrument would be the 
presentation to the operator and to the level control system 
of an indicated level inconsistent with the actual level. 



34 




Figure 5: Effective Reference Leg Height of Flooded 

Condensing Pot 



Condensing Pot 




Heat 
Rejected V 



States 

pressure p c 
average vapor volume 
fraction oi~ 



Total Volume = V c 
Total Mass = M c 
Total Enthalpy = H 



Figure 6: Condensing Pot Control Volume 

For the system under consideration, the steam generator, 
condensing pot, and the drain pipe, two control volumes are 
identified. The first consists of the condensing pot. The 

\ 
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second is the drain pipe itself (between the condensing pot 
and up to the steam generator tap connection) . They are 
depicted schematically in figures 6 and 7. 

The assumptions of the following calculations are: 

1. Each control volume is assumed to consist of saturated 
liquid and vapor in thermodynamic equilibrium under all 
operating conditions. 

2. The drain piping, usually consisting of various sections 
of horizontal and vertical piping runs, is simplified to be 
just one section of piping, inclined from the horizontal by 
an angle 0. This can vary from installation to installation. 

3. No significant concentration of non-condensible gases is 
present in the steam generator. 




Figure 7 : Drain Line Control Volume 



The model inputs are: p 9g (pressure in the steam generator 
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at the tap adjacent to the drain pipe connection) , p c 
(pressure in the condensing pot) , a c (the average vapor 
fraction in the condensing pot) , and h t (the liquid height 
in the drain pipe) . The model requires the initial values 
for p sg , p c , a c , and h t . The inclination of the drain pipe 
from the horizontal, 0, is assumed known. 

The key difficulty when trying to incorporate features of 
the drain pipe flow into the overall model is how to treat 
and calculate the transient liquid height in the tube when 
steam generator pressure changes. For the purposes of this 
model, the connection between the two is horizontal or 
nearly so and the liquid surface in the condensing pot is 
approximately a free surface. It therefore seems reasonable 
to conclude that the characteristics of the flow at the 
transition between the condensing pot and the drain pipe are 
best described by open channel flow theory. The liquid 
flow input is varied during transients in the sense that the 
liquid height in the pipe at the junction with the 
condensing pot is the result of weir flow. The liquid mass 
balance in the condensing pot is a result of weir flow and 
the instantaneous rate of condensation in the condensing 
pot. The liquid height in the drain pipe closer to the 
steam generator is a function of the liquid and vapor 
velocities and the pressure gradient along the pipe. 



4 . 2 Conservation Equations 

A first equation is obtained that defines conservation of 
mass in the condensing pot control volume. 



dM t 

dt 



C _ 



= Pviv^ 



( 4 ) 
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where 



M c 4 

Pv A 

jv A 

A 4 

Ih w 4 



mass of water in the condensing pot; 
density of water vapor in the condensing 
pot; 

superficial vapor velocity of the vapor 
entering the condensing pot from the 
steam generator; 

the cross sectional area of the drain 
pipe ; and 

the mass flow rate of liquid leaving the 
condensing pot in the weir. 



Note that the superficial liquid velocity in the weir is not 
necessarily the same as j ( , the superficial liquid velocity 
in the drain pipe close to the steam generator, at least 
during transient conditions when the pressure of the steam 
generator and the condensing pot are changing with time. 



Conservation of energy, in terms of enthalpy, in the 
condensing pot control volume gives : 



where H c 

V c 
A 

H v 

H, 

Q 

The total mass and enthalpy in the condensing pot are: 

M c =[p t (l-a c ) + p v a c ] V c (6) 



■^~Vc-J§-P v J v AH v -m w H r Q ( 5 ) 

4 the total mass-enthalpy product in the 
condensing pot; 

4 the total volume in the condensing pot; 

4 the total cross sectional area of the 
drain pipe; 

4 the enthalpy of the vapor entering the 

condensing pot from the steam generator; 

4 the enthalpy of the liquid leaving the 
condensing pot; and 

4 the heat loss from the condensing pot to 
ambient . 
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( 7 ) 



H c = [p t H t (l-a c ) +p v H v a c ] V c 

where p ( 4 density of the liquid leaving the 

condensing pot . 

In the drain line control volume, conservation of mass 
gives : 

-^=K + S> t J( A (8) 



where Mt 4 total mass of water in the drain line or 

tube and the influence of the vapor 
phase in the pipe has on the change in 
mass in the pipe has been neglected. 

As noted earlier, it is assumed that the transient flow 
condition at the entrance to the drain line, at the 
interface to between the condensing pot and the drain line, 
can best be modeled as a weir. From open channel flow 
theory, reference Cl and LI, for steep channels, the flow is 
critical at the channel inlet, thus the flow energy is 
minimized there, yielding the critical flow condition: 



3 tw 



=1 



(!-«*) 






g£> t cos© 



( 9 ) 



so that : 



lh w = p t j tw A=p t (l~ct b ) 



gD t c os0 



( 10 ) 



where 



jj« * 

Ob 4 

9 4 

g 4 



superficial liquid velocity at the weir, 
vapor fraction at the boundary between 
the condensing pot and the top of the 
drain pipe, 

angle of elevation of drain pipe from 
horizontal, 

acceleration of gravity, 
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D t 6 diameter of the drain pipe, and 

a« & energy or coriolis coefficient. 



Ob only applies in the vicinity of the weir and is not the 
same as the void fraction in the pipe, ot,., closer to the 
steam generator. The energy coefficient, a t , varies with 
the liquid depth and is usually determined empirically, is 
set equal to 2 for simplicity since, according to reference 
Cl, experimental values of a t are in the range of 1.03 to 
1.36 for fairly straight, prismatic channels. Since the 
precise degree to which the connection between the drain 
pipe and the condensing pot approximates an open channel is 
not known, the value for a ( is taken to be 2 as a worst case 
estimate . 

Differentiating (3) and (4) 




dM c 



dt 




and 



dH c JdH c \ 



dt dp dt [ da dt 
\ c / a e \ C 'P C 



dH e J d{p l H t ) {i _ a ) , d(p v H v ) 



dt dp c c dp c 
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substituting (11) and (12) into (4) and (5) 




dPc dP c 




+ (P v H v -p t H t ) V c -^=p v j v AH v -lh w H t -Q (14) 



The total mass in the drain line, M,., is: 

M t =A i L c p> t (15) 

where L t 4 the total length of the drain pipe; and 

A,= (l-a t )A. 

Relate Mt, the total mass of water in the drain line, to h { , 
the liquid height in the drain line, through the following 
geometric relationships (see figure 8) : 



i^OS ' 1 





(17) 



Substituting (16) into (15) : 




( 18 ) 
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Figure 8 : Relation Between h, and (]) in the Drain Pipe 



differentiating (14) 




L cPt 



dt 



- - cos 2 i 



dt 



+ sin 2 § t 




+ 2 sin 2 $ t 




- 2sin 2 § c 




(19) 
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recall 



-f-(cos' 1 x) = - — - — 

Using the above formula to differentiate (16) 

2 dh t 



t _ 



dt 



A 



“N 






dt 



( 21 ) 



substituting (21) into (20) 



dM t _D t 



dt 



L tPr 





( h, Y 




sin 2 


cos' 1 1-2 — - 


dh t 




1, M 



A 
N ^ 



i-A 

0. 



dt 



substituting (22) into (8) 



T' L e Pr 





( h, W 


sin 2 


COS' 1 ! 1 -2 —— 

[ a1 



A 
N ^ 



i- 






dh. 



Combining (13) , (14) , and (23) into a matrix 

d - 






(24) 



where : 



A= 



An -A.2 

■^21 -^22 



0 

0 



0 A 
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(25) 



and 



produces 



( 22 ) 



(23) 



equation : 
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Pc 




P vj v* ^ w 


x= 


«c 


( 26 ) £= 


P v jvAH v -lh w H t -Q 




h t 




K + Ptit A 



( 27 ) 



The entries of ^ are: 



*n = 



<1-«J + ' 



ap, 

dp c vx “ c/ "ap 

*12 — ( P v~ P < ^ ^ <. 






*21 = 






<P f ^i)|(l-a c ) + {a|- (PvH v )Ja c -l 

A 22 =(p v H v -p t H t )V c 



V. 



D, 

*33 -- 2 _i, cP« 



M , h 

N ^ c [ -DtJ 



4 . 3 The Functional Dependence of if 
4.3.1 i „ and i, 

While equation (27) identifies the components of £ used in 
equation (24) , it is useful for computation to clearly 
establish the functional dependence of the individual 
components of the 1 . The method for finding the 
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superficial velocities, j v and j c , is adapted from the 
approach used by Taitel and Dukler, reference Tl . Figure 9 
shows a smooth equilibrium stratified flow on which the 
following equations and geometric relationships are based. 

A momentum balance on the liquid and vapor phases yields : 




-A v -^-x wv S v -x 1 S i -A v p v gsinQ=0 (2 8) 

HP 

-A l ~-i wl S i + i i S i -A i p ( grsin0=O (29) 

where A,. k cross sectional area of the vapor phase; 

A* k cross sectional area of the liquid 
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wv m 

* 

ti A 

Sv A 

S ( A 

Si A 

0 A 



phase; 

wall to vapor phase shear stress; 
wall to liquid phase shear stress; 
interfacial shear stress; 
wetted perimeter of vapor phase; 
wetted perimeter of liquid phase; 
interfacial wetted perimeter; and 
angle of pipe elevation from horizontal 
(see figure) . 



using the nomenclature from T1 and C2 . 



From reference Tl, the geometric quantities are: 



sr 



where 



}=COS 



•f- 2 A 

l 



5 i =D c sin<J) 

S v =D c (:i-4>) 



A V =A-A ( 



The hydraulic diameters are: 



4 A v , 4 A. 

D v = - v and D t =—+ 

v s+s< * s v 



the velocities are: 



, JyA 

u. = — — and u=— — 



The liquid and vapor friction factors, using a Darcy 
Weisbach representation, are: 









and 




where 



Mv.i * 



dynamic viscosity of vapor and liquid, 
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n 



C. 



V, 



respectively; 

.184; and 

. 2 . 



for smooth stratified flow 




The shear stresses are evaluated from: 




_ 5PiMU| 



X 



wv 



_ v| ^v| 

8 



fiPvK—q|(Uy-q) 

8 



Equations (28) and (29) are 2 non-linear equations in 4 

. . dP 

unknowns: j v , , h, , and P' (where P* = -— ) . By defining j v 



equations (28) and (29) can be iteratively solved for j T and 
j t in terms of h, and P' . Note that P' is the difference 
between condensing pot and steam generator pressures divided 
by the length of the drain pipe. Equations (28) and (29) 
are iteratively solved for j„ and j { using known values of h 
and P' . Note, the equations used to solve for the 
superficial velocities are in error because the shear 
stresses above were calculated using a factor of 2 in the 
denominator, not 8. This error was discovered too late to 
re-run the simulations . The resulting sets of data were 
tabulated and used to correlate the superficial phase 
velocities given known values of liquid height in the drain 
pipe and pressure gradient along the pipe. The superificial 
velocity data is shown plotted in figures 10 and 11. 



and j ( as functions of h { and P' , that is: 

j v (P' ,h t ) and j t {P* ,h,) 
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Figure 9: Superficial Vapor Velocity Correlation 



Figure 10 : Superficial Vapor Velocity 
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Superficial Velocities 




Figure 10: Superficial Liquid Velocity Correlation 



Figure 11: Superficial Liquid Velocity 

4.3.2 The Weir Mass Flow Rate lh w 

ih w is known from equation (10) using the critical flow 

condition, (9) , but the relationship between a*, and a c , the 
vapor fractions at the condensing pot-drain line interface 
and in the condensing pot, respectively, must be more 
clearly identified. 

Note: (Xb is the vapor fraction at the boundary between the 

condensing pot and the top of the drain pipe, it is not the 
same as the vapor fraction throughout the drain pipe. 

The weir mass flow rate is: 

K=jiw A Pt do) 
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From the definition of a and figure 8, the following 
equations relate h t to 9: 



(1-a) 




(31) 



and 

A,=/-^| 2 (4>-cos<$>sin4>) (17) 



Note that these equations 
apply in the drain line and in 
the condensing pot 



Rearranging (16) 



2h >c , l 

— £=l-cos$ c 



(32) 



and solving for h tc (water height in the condensing pot) 

A <c =-| (l-cos<f> c ) (33) 



substituting (17) into (31) , and letting D=D C , <|)=<|) 0 , and a^a c 
yields 



| X)| 2 | <j> c -gi^4> e cog<}> c j 



(34) 



so that: 



A 

§ c -sin§ c costy c = (l-a c ) - — = (l-a c )it (35) 

( 4 ) 

If a c is known, (35) is used to find <|> c , then (33) is used 
to find h (o . Knowing h tc , h w (liquid height at the weir, see 
figure 12) is found by: 
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If 



then h w = 0 (36) 




Figure 12: Relation Between h, e and 

Where the small angular dependence of h w is neglected. 
Clearly, only a certain range of h {c values will produce 
values of h w in the range from 0 to D t . If h tc is above the 
top of the drain pipe, the pipe entrance will always be 
filled. If h (o is below the bottom of the drain pipe 
connection to the condensing pot, the pipe will be empty of 
liquid. 

Once h w is found, 0^, the vapor volume fraction at the weir, 
is computed from: 

Now that is known, ih w , the liquid mass flow rate at the 
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(39) 



(4Q) 

A t 11 



weir present as the condensing pot-drain pipe interface, is 
found using equation (10) : 



(I - ®*) 




( 10 ) 



4.3.3 The Heat Transfer Rate Q 

The heat lost from the condensing pot to the reactor 
building can be expressed as : 

Q=UA(T gCm -T aiz ) (41) 



UA is calculable from the geometry of the condensing pot and 
some empirical relations for the film heat transfer 
coefficients of the liquid condensing on the inner surface 
of the condensing pot and of the air on the outer surface. 



For a cylindrically shaped condensing pot. 




± \ + ± 

izr^c^stm 2nL c ic g£; 2 TZI out L c h a ^ I 



where r in 4 inner radius of the condensing pot; 

r out 4 outer radius of condensing pot; 
h 3tm 4 film heat transfer coefficient for the 
water vapor condensing on the inside 
surface of the condensing pot; 
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A 



A 

A 



film heat transfer coefficient for air 
in the reactor building; 

thermal conductivity of condensing pot; 
and 

inner length of condensing pot . 




Figure 13 illustrates the typical temperature profile from 
the condensing pot to ambient. An empirical expression for 
h etm , from reference M3 based on the Nusselt method for 
predicting the rate of heat transfer from a vapor condensing 
in a laminar film and running down a vertical surface, is 
used to get an order of magnitude estimate for the size of 
the film heat transfer coefficient of the condensing water 
vapor : 






where Nu, 



J h fg p t gL c 
' {^k(T atm -T la ) 



(43) 



evaporation increment; 
thermal conductivity of water; 
dynamic viscosity of water; 



where : 



k 



A 

A 

A 
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L, 



C 



A 



g 

p 



a 

a 



T atm A 



A 



wetted length of condensing pot above 

normal water level; 

acceleration of gravity; 

density of liquid water at condensing 

pot pressure; 

Temperature corresponding to saturation 
pressure in condensing pot; and 
Temperature of inner surface of 
condensing pot; 



assuming : 



1. Vapor which flows parallel to the surface flows so slowly 
that its drag on the condensed film can be neglected. 

2. Motion of any liquid element in the condensed film is 
governed solely by a balance of the gravity and viscous 
forces acting on the element, hydrostatic pressure forces 
and inertia forces can be neglected. 



Despite the fact that the Nusselt correlation applies 
specifically to vertical plates, it is useful for getting an 
estimate of the magnitude of h 3tm . 



To estimate h air , a simplified expression for the heat 
transfer coefficient is used, obtained from reference Ml, of 
the form: 

<44) 

where A and b are constants, depending on geometry and flow 
conditions, and L, is the significant length, also depending 
on type of geometry and flow. Using the values suggested by 
reference Ml for an horizontal cylinder and laminar flow, 
A=.27, b=l/3, L a is the diameter of the condensing pot, and 
At is T 0 , the outer surface temperature of the condensing 
pot, minus T^, the temperature of the ambient air in the 
reactor building. Heat transfer by radiation from the 
condensing pot to the reactor building, a factor that is 
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less than 20% of the overall heat transfer coefficient of 
air, is ignored. Including the contribution due to 
radiative heat transfer would slightly increase the value of 
the calculated h air . 



The difficulty in computing both the film heat transfer 
coefficients is that the value of each coefficient depends 
on the magnitude of the temperature difference across the 
respective films, which one cannot know until the heat 
transfer coefficient is known. The heat transfer rate from 
the condensing pot to the air of the reactor building, 
assuming that an estimate of j t is available to compute ihj , 
can be found from: 

Q=m t h fg (45) 



where: ih t A mass flow rate of liquid condensing in 

the condensate pot . 



If the thermal resistance of the condensing steam film is 
known, the temperature difference across the film can be 
expressed by: 






_ Q 



R ♦ 



and R 



2 nr ln h gCm L c 



where R r A thermal resistance of the condensing steam 

L sca 

film, which is a function of the temperature 
difference across the film. 

Equation (46) can be solved iteratively to find the thermal 
resistance of the condensing steam film and thus h stm . The 
same procedure is used to find h air . 

Figure 14 shows a typical condensing pot and drain pipe 
arrangement. Using a pressure of 7.34 MPa to compute h £g , 



55 



and an estimate of .02 for j t , ih 2 and Q can be computed 

using equations (10) and (45), respectively. The wetted 
length, L c , in equation (43) is taken as one half of 3.5 



inches, the length of the condensing pot shown in figure 14, 
since the condensing pot is nominally half filled with 
liquid at all times. 




Figure 14: Geometry Of Condensing Pot 

Taking k = .626 W/ (m-K) , g = 9.807 m/sec 2 , L c = 44.5 mm, r in 
= 24.6 mm, (1=90 [iPa-s, p x = 729 kg/m 3 , an iterative 
solution of equation (46) and using equation (43) for h stm 
yields a value of h 9tm = 5400 W/ (m 2 -K) . An iterative 
solution of equation (46) using the expression for h air of 
equation (44) and assuming that h stm is so large in 
comparison with h air that its effect of the total 
differential temperature is negligible, gives h air = 12 
W/ (m 2 -K) . Returning to equation (42) , using the above 
values and r out =27.4 mm, k st = 50 W/ (m-K) , one finds: 
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1 



=0 . 000663 



^ Z ’ ln^~ 'c^ 'steam 



In- 



■ in 



■ out 



2 u L c k gc 



=0 . 008 



^ out^c^a 



=10.09 



These results clearly show that the air film heat transfer 
coefficient dominates in the denominator of equation (42) . 
For the purposes of simulation, a constant value of h air is 
used because the actual value during transients at normal 
operating pressure, a function of the pressure in the 
condensing pot, will not vary too far from the nominal value 
calculated above . 



The resulting functional dependence of Q is 



Q=UA(T sCm (p c )-T air ) (47) 



Having specified the functional forms of each component of 
? from equation (24) , it is easier to see how to compute 
the necessary partial derivatives that will be needed to 
solve (24) . 



4 . 4 Calculations 



The equations to be solved are of the form: 




1 



and 
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f= 



fli 

f2 

f3 



where : 



f^pyiPc) 3 V {P' ,h t ) A-m w {a c ) 



f 2 = p v (p c ) j v (P', h t ) H v {p c ) A-m w (a c ) H t (p c ) -£>(p c ) 
f 3 =ih w {a c ) +p, (p c ) j ( {P', h ( ) A 



and 



P ‘ 



= Jp £ -p ££ ) 

Lt 



Converting to a one-step, implicit difference of the 
linearized equations: 



x n ^-x n 



At 






||j (x n+1 -x n ) 



and we seek x 11 * 1 : 



x a * 1 = 



A n _ 


d?) D 


At 


[dZ) J 



i-i 



£ D +x n 



(48) 



4.4.1 The Jacobian of the f vector 

The Jacobian of 1 is computed using numerical 
differentiation, where the Jacobian component in row i 
column j is: 
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( aif'l _ fjpCj+AXj) -fj(Xj) 
axj 



Of course, it is important to choose the A value small 
enough to avoid errors due to nonlinearity of the subject 
functions and large enough to avoid roundoff in getting a 
meaningful average value for the derivative. Several plots 
of the functions f lt f 2 , and f 3 were computed to estimate 
adequate A values. The values used are: 

4P c =30 Aa c =.0005 Ah j = .0005 



The complete Jacobian, 



d? 

dx 



of 1 is : 







3*x 


3*x 


9f il 


d Pc 


da c 


dh 1 


df 2 


3*2 


3*2 


dPc 


da 

c 


dh 1 


df 3 


Mi 


3*2 


dp c 


da c 


dh 2 



For details on the software used to solve equation (41) for 
various initial conditions of x and pressure transients in 
the steam generator, please refer to the Appendix. 

4.4.2 Geometry and Initial Conditions 

For the purposes of these calculations, the following 
average parameters were used: 

Length of drain piping (L t ) =2.5 m; 

Reactor building ambient air film heat transfer coefficient 
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(h air ) =11 W/(m 2 -K); 

Inner radius of condensing pot (r in )=24.6 mm; 

Outer radius of condensing pot (r out ) =27.4 mm; 

Ambient temperature of Reactor Building (T. mh ) =50 °C 

To calculate the initial or steady state values of x, a 
initial value for p c is assumed. Knowing this variable and 
the ambient reactor building temperature, equation (46) is 
used to find Q, the heat rejected by the condensing pot to 
ambient. With Q known, the liquid mass flow rate from the 
condensing pot the steam generator is found from equation 
(45) . 

Equation (10) is used to compute the superficial phase 
velocities, presuming that the liquid and vapor mass flow 
rates are equal at steady state. Once j v and j ( are known, 
equations (28) and (29) are solved iteratively until values 
for h and P' are found. With a calculated P' and a known 
value of p sg , p c can be computed and the process repeated 
until convergence is obtained for a value of p c . This is 
how the steady state values of x, composed of p c , a c and h* 
were computed. This process is diagrammed schematically 
below: 



The above procedure for computing the initial conditions 
based on a combination of steam generator pressure, p sg , an 
assumed pressure gradient between the condensing pot and 
steam generator, P' , and the geometry of the condensing pot. 



P C Q b w ® c 
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It is, however, possible to calculate a liquid mass flow 
rate from the condensing pot to the steam generator, m w , 

that makes it impossible to solve the momentum balance 
equations for the liquid and vapor phases in the drain pipe, 
equations (28) and (29) . Recall that at steady state, 
liquid superficial velocity in the drain pipe, j lr is 
related to ifl w strictly by the liquid density and the pipe 

cross sectional area (equation (10) ) . Reference W1 notes 
that a plot of j t and j v for constant liquid height, h { , 
figure 23 for the tabular data calculated, enables one to 
identify graphically the regions of permissible operation in 
the drain pipe. Positive and positive represents 
cocurrent flow up from the steam generator to the condensing 
pot using the positive sense shown in figure 7. Negative j„ 
and negative j t represents cocurrent down flow from the 
condensing pot to the steam generator. The region of 
positive vapor flow and negative liquid flow between the 
axes and the locus of points joining tangent lines to the 
lines of constant vapor liquid height, called the flooding 
line, represents the limited region of counter current flow, 
where solutions are possible (see figure 23) . Operation 
above and to the left of the flooding line is not possible, 
because the drain pipe is already flooded at that point. 

If, while calculating the steady state solution, one picks 
assumed conditions of pressure gradient, P' , that results in 
superficial velocities that are beyond the flooding line, 
one must continue to select a lower pressure gradient and 
recalculate the superficial velocities until physically 
realizable values of the superficial velocities are 
obtained. 



4.4.3 Model Limitations 
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The tabular data, values of j„ and j t for known values of h, 
and P' , must be calculated before implementing the other 
model equations. The limits of the allowed liquid heights 
in the drain pipe are determined by the drain pipe diameter, 
but the range of potential pressure gradients cannot be 
known until after the model is implemented and some 
solutions calculated. The possibility exists for steam 
generator pressure changes to occur rapidly enough for the 
pressure gradient between condensing pot and steam generator 
encountered during transients to exceed the maximum or 
minimum values of P' used to compute the tables of j { and 
data. To account for this possibility, there are two 
options. First, perform extensive tests of the model using 
the worst case positive and negative steam generator 
pressure transients and expand the range of drain pipe 
pressure gradients used to generate j v and j, by trial and 
error. This is fairly easily, though tediously, 
accomplished and only increases the memory that the smart 
instrument must have. Second, if the smart instrument is 
capable of performing the calculations rapidly enough, the 
designer could elect to solve equations (28) and (29) at 
each time step rather than relying on tabular data computed 
in advance . 

4.4.4 Solutions 

The independent variable in the implicit difference equation 
(45) is p 9g , steam generator pressure. Equation (48) is 
solved for the state variables condensing pot pressure, p c , 
condensing pot vapor volume fraction, a c , and drain pipe 
liquid level, h t , for both constant steam generator pressure 
and up and down steam generator pressure ramps. 

4. 4. 4.1 Constant Pressure 
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Table 1 show the results of the solutions computed, at a 
steam generator pressure of 7.34 MPa and various condensing 
pot and drain pipe diameters . 
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Steam Generator 
Pressure = 7.34 


MPa 




Pet 

Chg 




Pet 

Chg 


Condensing Pot 
Diameter (mm) 


50.8 


50.8 


0 


54 . 8 


7 . 9 


Drain Pipe 
Diameter (mm) 


25.4 


27 . 4 


7.9 


27 . 4 


7 . 9 


(P=-P sg ) (Pa) 


-183.19 


-143.16 


21.9 


-143.16 


21.9 


Condensing Pot 
Vapor Volume 
Fraction a c 


.7123 


.741 


4.0 


.724 


2.3 


Drain Pipe Liquid 
Level h { (mm) 


5.47 


4 . 84 


11.5 


4 . 84 


11.5 



Table 1 : Steady State Solutions for Constant Steam 

Generator Pressure 

A nominal steam generator pressure of 7.34 MPa was chosen 
since many commercial steam generators operate between 7 . 6 
MPa and 6.9 MPa when steaming at normal operating pressure. 
The piping dimensions are varied to see the influence of 
dimensional tolerances on the stability of equation (48) . 

As the data in table 1 shows, P' is the most sensitive state 
variable with respect to the diameter of the drain pipe. A 
7.9 percent increase in drain pipe diameter reduced the 
pressure drop by 21.9 percent, the vapor volume fraction 
rose by only 4 percent and the liquid level decreased by 
only 11.5 percent. Physically, it makes sense that 
enlarging the drain pipe diameter would reduce the wall 
friction forces, thus reducing the pressure drop. The vapor 
volume fraction would also decrease since expanding the pipe 
diameter reduces the liquid level in the pipe. The last 
column of table 1 indicates that condensing pot vapor volume 
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fraction is the only state variable sensitive to condensing 
pot diameter and only weakly at that. One should also note 
that the pressure drop between the condensing pot and the 
steam generator is an extremely small fraction of steam 
generator pressure, about .0025 percent. 

4. 4. 4. 2. Pressure Ramps 

Table 2 shows the steady state solutions reached for equal 
and opposite steam generator pressure ramps of 300 kPa, both 
starting from 7.34 MPa. Pressure ramps were used as a first 
approximation for the generator pressure changes that result 
from changes in turbine loading. 







Down Ramp 


l 

Up Ramp | 

i 


Pressure 

• 

Ramp (kPa) 


0 


-300 


[ 

+300 

! 

1 


(Pc-Pag) (Pa) 


-184.19 


-184.46 


-179.42 


Condensing 
Pot Vapor 
Volume 
Fraction a c 


.7123 


.7144 


.7102 


Drain Pipe 
Liquid Level 
h { (mm) 


5 . 47 


5.44 


5.51 



Table 2 : Steady State Solutions for Pressure Ramps 

The trend of the steady state data corresponds well to what 
is predicted by the physical situation. For the negative 
pressure ramp, lower pressure in the steam generator results 
in less condensation in the condensing pot and reduced heat 
transfer from the condensing pot to ambient since the 
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condensing pot is assumed to be saturated and thus the 
temperature of the water vapor in the condensing pot is 
lower for the lower pressure. Reduced condensation in the 
condensing pot results in a lower liquid level there which 
causes the liquid level in the drain pipe to be lower (from 
equations (36) to (38) ) . The opposite trend is observed for 
the positive pressure ramp. 

As mentioned in section 4.5.3, since the values of P' used 
to solve equations (28) and (29) for and j { were 
specified in advance of computing the solution to equation 
(48) , the possibility exits that transient values of P' 
outside the ±600 Pa/m range those used might result. 
Preventing this form occurring constrains the maximum 
allowable rates of steam generator pressure change that my 
model could accept. As shown in table 3, the maximum 
allowed steam generator pressure gradient with respect to 
time is also different for up and down pressure transients . 



Generator Pressure Ramp 


Maximum Rate of Change 


Down 


1.33 kPa/sec 


Up 


.42 kPa/sec 



Table 3 : Maximum Generator Pressure Gradients 

For the negative pressure gradient, figures 15 through 18 
show the transient steam generator pressure, the 
differential pressure between the condensing pot and steam 
generator, the vapor volume fraction, and liquid height in 
the drain pipe. Note that both the differential pressure 
and the drain pipe liquid height asymptotically approach 
constant values during the transient, but rapidly converge 
toward their steady state values after the transient. It is 
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particularly interesting to note that the these solutions to 
equation ( 41 ) show that the steam generator level model is 
capable of directly inferring the proximity to a flooding 
condition because an instantaneous value of liquid height is 
calculated during the transient. Other methods of 
predicting flooding rely on correlations between the 
superficial liquid and vapor velocities in the pipe since 
liquid height in the pipe is not usually calculated. 



xlO 6 -300 kPa Pressure Change In 3.75 Minutes 




Figure 15: -300 kPa Pressure Ramp, p, g 
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Differential Pressure for -300 kPa Pressure Change In 3.75 Minutes 




Figure 16: -300 kPa Pressure Ramp, p e -p, g 

xlO* 3 Vapor Volume Fraction for -300 kPa Pssr Chg In 3.75 Min 




Figure 17 : -300 kPa Pressure Ramp, a c 
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xlO ' 3 



Drain Line Liquid Level for -300 kPa Pssr Cluj In 3.75 Min 




Figure 18: -300 kPa Pressure Ramp, h, 

Figures 19 through 22 show the transient steam generator 
pressure and state variables for the positive steam 
generator pressure ramp. As discussed, the positive 
pressure ramp had to occur over a time interval 3 times 
longer than that for the negative ramp in order to reach a 
solution. It is interesting to note that for the positive 
pressure transient, the transient liquid height in the tube 
increases just as it did for the negative pressure 
transient . 
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Psg (Pa) 



xlO 6 +300 kPa Pressure Change In 13.33 Minutes 




Figure 19: +300 kPa Pressure Ramp, p, g 
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Differential Pressure for +300 kPa Pressure Change In 13.33 Minutes 




Figure 20: +300 kPa Pressure Ramp, p c -p, g 

xlO 3 Vapor Volume Fraction for +300 kPa Pssr Chng In 13.33 Min 




Figure 21: +300 kPa Pressure Ramp r a c 
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h J (m) 




Figure 22 : +300 kPa Pressure Ramp, h, 



For the negative pressure ramp, the rising trend in drain 
pipe liquid height can be explained as follows: 

1. (p c p s g) changes from negative to positive as 
soon as steam generator pressure begins to decrease 
because condensing pot pressure cannot change as 
rapidly as steam generator pressure. 

2. To balance the momentum equations, (28) and (29), 
j v must become less than zero (see figure 10, j v must 
be less than zero when p c -p 3g becomes greater than 
zero) . 

3. Since the initial change in p c is small, a c remains 
constant or changes slowly at the start of the 
transient, suggesting that h { is initially constant. 

4. Looking at figure 23, if j v becomes negative and h ( 
remains constant, must become more negative as well, 
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possibly by a factor of 2 to 3 . 

5. Figure 11 shows that as j x becomes more negative 
for a fixed value of P' , h x rises. 

For the positive pressure ramp, the rising trend in drain 
pipe liquid height can be explained as follows: 

1 . p c _ P a g becomes even more negative as soon as steam 
generator pressure begins to increase. 

2. To balance the momentum equations, (28) and (29), 
j„ must become more positive (see figure 10, must 
increase when p c -p sg becomes more negative) . 

3. Since the initial change in p c is small, a c remains 
constant or changes slowly at the start of the 
transient, suggesting that h { is initially constant. 

4. Looking at figure 23, if j v becomes more positive 
and h x remains constant, must become less negative. 

5. Figure 11 shows that as j x becomes less negative 
for a fixed value of P' , h ( rises. h x would actually 
decrease if the P' only became slightly more negative. 
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One Dimensional Horizontal Flow-const h 




Figure 23 : Flow Map for Drain Pipe 



4 . 5 Flow Transitions 

A major goal of the steam generator water level model is to 
create an instrument sophisticated enough to alert the 
operator of the approach to flooding conditions in the drain 
pipe. If flooding does not occur, the pressure difference 
across the drain pipe remains very small so no correction is 
required. Recall that the drain pipe returns the overflow 
of condensed steam from the condensing pot to the steam 
generator and should be sized large enough to transport the 
maximum expected condensation flow without undergoing a flow 
transition from stratified flow to annular or intermittent 
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flow. Such a flow transition changes the height of the 
re f® r snce leg and produces erroneous level readings . 

4.5.1 Taitel and Dukler 

Taitel and Dukler (reference T2) developed a theoretical 
model for predicting the onset of flow regime transitions in 
near-horizontal piping that agrees well with experimentally 
observed results. 



The procedure was first to visualize a stratified liquid and 
then the mechanism describing the change from stratified 
flow. Beginning with stratified flow, as the liquid 
velocity increases, the liquid level rises and a wave is 
formed which grows rapidly, tending to block the flow. At 
lower vapor velocities, the wave crest forms a complete 
bridge and slug or plug flow results . At higher vapor 
velocities, there is insufficient liquid flowing to maintain 
or even form the bridge and the liquid in the wave is swept 
up around the pipe to form an annulus with some entrainment 
if the vapor velocity is high enough. 



Taitel and Dukler found that the criterion for the 
transition from stratified to intermittent and to annular 
dispersed flow is: 



dA t 




(48) 



where F is the Froude number modified by the density ratio: 



F=. — ^ J v . - - (49) 

\ Pi~ Pv ^/D^gcosB 



C2 is estimated as l-h ( /D t and the dimensionless values in 
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equation (48), signified by a ~, are found as follows. 






fir 



h 



A v «-J[c OS" 1 ( 2£ t ) - ( 25,-1 ) ft - (2fi,-l ) 2 ] 
COS' 1 (2£ { ) - (2^-1) yjl- (25,-1) 2 ] 



.A— 




Clearly, if one knows the liquid height in the drain pipe 
and the superficial vapor velocity, j„, the proximity to the 
Froude number characteristic of the flow transition can be 
calculated from equations (48) and (49) . The chief obstacle 
to directly implementing a scheme to evaluate the Froude 
number is that the superficial velocities and liquid height 
in the drain pipe between the condensing pot and the steam 
generator are not measured. Despite the numerical error 
that was made when computing the shear stresses in equations 
(28) and (29), the results of section 4.4 suggest that it is 
theoretically possible to infer these variables from 
knowledge of the steam generator pressure and piping 
geometry . 
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4.5.2 Drift Flux Model 



As discussed, the flow in the drain pipe is normally two 
phase, counter current flow, with steam flowing from the 
steam space of the steam generator to the condensing pot and 
condensate returning from the condensing pot to the steam 
generator, illustrated in figure 9. The drift flux model, 
presented in reference Wl, a separated flow model focusing 
on the relative motions of the individual phases in the 
pipe, is well suited to analysis of bubbly and slug regimes 
of gas-liquid flows and is appropriate to the situation at 
hand. 

The limit of counter-current flow, called flooding as the 
pipe fills completely with liquid, predicted by the drift 
flux model can be seen on a plot of solutions of the 
momentum balance equations, (28) and (29) , for constant 
liquid height in the drain pipe. Such a curve is plotted in 
figure 23, with superficial liquid velocities plotted 
horizontally and superficial vapor velocities plotted 
vertically, with each curve representing a different drain 
pipe liquid height. The region of negative liquid velocity 
and positive vapor velocity, called the limited region of 
counter current flow in reference Wl, is bounded by the axes 
and the dotted line known as the flooding line. The 
flooding line is the locus of points tangent to the j { , j v 
curves in the limited counter-current flow region, since an 
increase in the magnitude of either phase velocity beyond 
the flooding line is condition for which no steady flow 
solution is possible and a flow transition must occur. 

Conclusions 



The suggested model, reflecting the influence of important 
physical features for processes influencing the measurement, 
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produces solutions that are stable and finite for both 
constant steam generator pressure and steam generator 
pressure gradients. Although a numerical error was made 
when solving the two phase flow equations (equations (28) 
and (29) ) , the steady state solutions appear to correspond 
well to the physical situation. The solutions produced for 
steam generator pressure ramps indicate that a potential 
shortcoming of the model is that some transients might cause 
the differential pressure between the condensing pot and the 
steam generator to exceed the values of differential 
pressure used to pre-compute superficial liquid and vapor 
velocities in the drain pipe. 
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Chapter 5 



Smart Instrument Applied to Steam 
Generator Level 



Having postulated a steam generator level measurement model 
capable of determining the pressure in the condensing pot 
and the liquid level in the drain pipe as functions of time, 
is appropriate to consider how successful a smart instrument 
including such a model might be at improving the accuracy of 
level measurement over a standard level sensing device. 

Since no operational plant level data is available for 
comparison, the Prism (reference Kl) program is used to 
generate steam plant data for an up power and a down power 
transient at high power. 

5 . 1 Sensor Dynamics, Information 

The relationship between steam generator water level and the 
differential pressure cell output is: 



Equation (51) is the same as equation (2) , except that the 
functional dependence of all the variables has been 
specified. Thus the liquid density, p dc , is a function of 
the temperature in the downcomer, the vapor density, p v , is 
a function of the steam generator pressure, and the 
reference leg liquid density, p r , is a function of the 
ambient temperature in the reactor building. 




(51) 
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For a conventional level instrument, we can make some 
simplifications to make the level easier to compute. Over 
the pressure range usually encountered in commercial power 
reactors between 0 and 100 percent power, the vapor density 
changes from 36.5 kg/m 3 to 39.5 kg/m 3 so using an average 
value of 38.0 kg/m 3 for the vapor density introduces an 
error of no more than 4 percent . Assuming that the reactor 
building ambient temperature does not change, the reference 
leg density can be regarded as constant at 988 kg/m 3 . The 
conventional instrument can also neglect the pressure 
differential between the generator and the condensing pot. 

5 . 2 Downcomer Temperature Effects 

As noted in the description of steam generator internals, 
the steam generator downcomer is the region between tube 
sheet wrapper and the generator vessel, below the feedwater 
inlet ring. Prior to entering the tube bundle area, 
feedwater returning from the condenser mixes with saturated 
liquid that has been trapped by the moisture separators in 
the downcomer. The saturated liquid from the moisture 
separators also preheats the feedwater. The density of the 
downcomer liquid is related to its temperature. 
Unfortunately, downcomer temperature is not a measured plant 
variable, although feedwater temperature is. To estimate 
downcomer temperature, a heat balance is performed on the 
downcomer : 



lh dc C p T dc = lh fu C p T fw + (lh dc -lti fw ) C p T gaC (52) 



where 



Aic a 

c p A 

T dc A 

ih fw 4 

T fw a 



mass flow rate in rate in downcomer; 
specific heat capacity of water; 
temperature of downcomer water; 
mass flow rate of feedwater; 
temperature of feedwater; and 
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T aaC a saturation temperature of generator 
pressure . 

assuming that the specific heat capacities in equation (52) 

i 

are approximately equal. Rearranging (52): 

rr _ + dc ^ dc ^ ^ sat / C \ 

T dc * (53) 



The feedwater and steam mass flow rates of the steam 
generator are measured plant variables, but the downcomer 
mass flow rate is not . The downcomer mass flow rate can be 
estimated using the Circulation Ratio. The Circulation 
Ration is a measure of the quantity of water entering the 
tube bundle section of the steam generator (downcomer flow) 
compared tot he total feedwater flow. Circulation ratio is 
defined as: 

CR _ Mass Flow Rate Entering the Tube Bundle _ ™dc 
Mass Flow Rate of Feedwater m fw 



Typically, the Circulation Ratio decreases from about 25 at 
lower powers to about 5 for high powers. By estimating the 
Circulation Ratio for the power level and multiplying it by 
the feedwater mass flow rate, it is possible to generate a 
reasonable value for downcomer mass flow rate. 

For the conventional instrument, it is assumed that a good 
estimate of downcomer temperature would be made and used to 
calculate a constant downcomer density. It turns out that 
the downcomer temperature is not a very sensitive function 
of feedwater temperature. For a feedwater flow rate of 20 
kg/s, feedwater temperature of 132 °C, saturation temperature 
of 290 °C, and a Circulation Ratio of 20, (corresponding to 
about 5 percent power), the downcomer temperature is 282°C. 
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For a feedwater flow rate of 400 kg/s, feedwater temperature 
of 268°C, saturation temperature of 286°C, and a Circulation 
Ratio of 5 (corresponding to nearly 100 percent power) , the 
downcomer temperature is 282 °C. Clearly, the conventional 
instrument is capable of using a single value for downcomer 
liquid density, 747 kg/m 3 based on T dc =282°C. 

A smart level instrument could accept feedwater flow rate 
and feedwater temperature as inputs and, with knowledge of 
the details of Circulation Ratio as a function of power, 
calculate an instantaneous downcomer temperature using 
equation (53) . This calculation is performed by the smart 
instrument in this report. 



5 . 3 Inputs and Outputs 

Figures 24 and 25 show block diagrams for the conventional 
and smart instruments. The conventional instrument's only 
variable input would be the output from the differential 
pressure cell. The only output for the conventional 
instrument is level, L d . The smart instrument would have 
the differential pressure cell input as well as the 
temperature of the feedwater, T fw , the mass flow rate of the 
feedwater, ih fw , the Circulation Ratio appropriate for the 

power level, CR, the steam generator pressure, p 3g , and the 
temperature of the reactor building, T^. The smart 
instrument would provide two outputs: level, L, and a 

flooding warning to alert the operator of an approach to 
flooding in the drain pipe . 
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Figure 24: Conventional Level Instrument Block Diagram 
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Figure 25 : Smart Level Instrument Block Diagram 
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5 . 4 Calibration 



The output data from the Prism program after running 30 
minute simulations of up and down power transients is used 
to calibrate the conventional and smart instruments . For 
the purposes of comparison, this procedure assumes that the 
steam generator data recorded by Prism represents the actual 
value or truth. The calibration procedure was to compute 
the downcomer water temperature as a function of time using 
equation (55) , assuming a Circulation Ratio of 5 for the 
power level, to find the downcomer liquid density, p t , as a 
function of time, compute the vapor density, p v , as a 
function of time based on the Prism data for generator 
pressure, use fixed values of reference leg density, p r =988 
kg/m 3 , reference leg height, H=2 m, and use the Prism data 
for generator level to solve equation (1) for the 
differential pressure cell output, aP, as a function of 
time . 

Naturally, using this calculation procedure, it is not 
instructive to compare the smart level instrument calculated 
level output with the level signal recorded by Prism. They 
are identical because the procedure and data used by the 
smart instrument to calculate steam generator level from 
differential pressure cell data is just the inverse of the 
procedure used to derive the differential pressure cell data 
from the Prism level information. 

5 . 5 Comparison 

The time varying differential pressure cell output, aP (t) , 
extracted from Prism data by the procedure described in the 
previous section, was applied to the simplified calculation 
(51) performed by the conventional instrument and to the 
more detailed calculation performed by the smart instrument. 
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Comparison for Down Power 




Figure 26: Comparison of Level Signals for Down Power 

Figure 26 is plot of the level signals as functions of time 
for the smart and conventional instruments for a down power 
from 100 to 90 percent power. Recall that the difference 
between the two calculated levels is that the conventional 
instrument uses an average value for vapor density, constant 
downcomer liquid density, and neglects the differential 
pressure between the condensing pot and steam generator. 

The smart instrument calculates the vapor density and liquid 
density at each time step, and does not neglect the pressure 
difference between steam generator and condensing pot. 

Figure 27 show that the level error between the conventional 
and smart instruments varies, but is generally below .5 
percent . 
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Level Error - Down Power 




Figure 27 : Level Error For Down Power 



Figure 28 is a comparison of the level signals for an up 
power transient from 90 to 100 percent. From figure 29, the 
minimum deviation between the level signals for the up power 
is .4 percent, but at no time is the error greater than 1 
percent . 
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Figure 29: Level Error For Up Power 



Conclusions 



The operation and accuracy of a smart instrument was 
investigated. The model was used by the smart instrument 
was developed in chapter 4 . The input variables for the 
smart instrument are steam generator pressure, feedwater 
flow rate and temperature, reactor building temperature. 
Circulation Ratio, and differential pressure cell output. 

For the basis of comparison, a "conventional" level 
instrument was developed by analyzing the behavior of the 
variables in equation (51) over the range of steam generator 
pressures reasonable during power operation and taking all 
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quantities in equation ( 51 ) constant except for the 
differential pressure cell output. The only input variable 
to the conventional instrument is thus differential pressure 
cell output. Since not operational differential pressure 
cell data was available, the Prism program was used to 
generate level information as a function of time for two 
power transients and the differential pressure cell data for 
both transients was extracted from the Prism level data. To 
the extent that the simulated differential pressure signal 
is valid, the smart instrument was more accurate than the 
conventional instrument, but the difference between the two 
level signals over the thirty minutes of simulation never 
exceeded one percent. It is likely that process noise, 
normal sensor fluctuations, and meter error would completely 
obscure such a small accuracy improvement. 

Although the smart instrument is not significantly more 
accurate than the postulated conventional level instrument 
over the range of pressure changes produced by the small 
power transients (ten percent) imposed on the Prism program, 
the smart level instrument is capable of estimating the 
liquid level in the drain pipe between condensing pot and 
steam generator while the conventional instrument is not. 
Thus, the smart instrument, while being only slightly more 
accurate, is capable of identifying and warning the operator 
of a violation of the assumption common to both smart and 
conventional level instruments : that the reference leg 

height is constant. While a utility's response to a smart 
instrument reporting that drain pipe flooding during power 
operations would no doubt be exactly the same as if a 
conventional level instrument's output produced unexpected 
jumps in level characteristic of drain pipe flooding (i.e., 
replace the drain pipe or condensing pot with different size 
piping), the smart instrument's ability to infer variables 
not easily measured and take action based on that 
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information has been demonstrated. 
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Chapter 6 



Conclusions 



A procedure for designing and implementing a smart 
instrument was applied to a Pressurized Water Reactor steam 
generator water level measurement. Important physical 
effects on generator level were considered and applied to 
find an overall model capable of inferring desired variables 
that are not directly measured. The model is stable for 
constant steam generator pressure and is also stable for a 
limited range of steam generator pressure transients . A 
conventional level instrument, using a simplified procedure 
to calculate level, was chosen for comparison with the smart 
level instrument. Both instruments received differential 
pressure cell data extracted from Prism level for two 
transients. As expected for this comparison, the smart 
instrument output conformed more closely to the "true" steam 
generator level, but the difference between the smart 
instrument and the conventional instrument was never more 
than 1 percent . 

6 . 1 Significance of Corrections 

To understand why there was so little difference between the 
smart instrument level data and the conventional instrument 
data, it is necessary to compare the magnitude of the 
corrections made by the smart instrument to the differential 
pressure signal. 

The differential pressure cell output extracted from Prism 
output data is on the order of 365 kPa. The maximum steam 
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generator to condensing pot pressure difference calculated 
was 1 kPa. Clearly, it is safe to conclude that the 
differential pressure is negligible. 

The smart instrument calculates the downcomer temperature 
and uses this temperature to adjust the liquid density in 
the downcomer at each time step. The conventional 
instrument uses a constant downcomer temperature of 282°C. 
The range of downcomer temperatures calculated by the smart 
instrument is 280°C (p t =750 kg/m 3 ) to 282°C (p t =747 kg/m 3 ) . 
This small variation in downcomer temperature produces a 
variation in liquid density of only 0.49 percent. 

The steam generator pressure changes for both the up and 
down powers were on the order of 100 kPa from a nominal 
steam generator pressure of 7 MPa. This pressure deviation 
results in a maximum vapor density difference over the 
entire maneuver of 1.6 percent. 

Originally, the temperature of the reactor building was 
intended to be an input to the smart instrument since it 
seemed plausible that a steam leak in the reactor building 
would significantly raise the ambient temperature in the 
vicinity of the reference leg, thus decreasing the density 
of the liquid in the reference leg. Indeed, a steam leak 
that raised the reactor building temperature from 50 °C to 
200°C would change the density of the reference leg from 988 
kg/m 3 to 8 65 kg/m 3 , a change of 12.5 percent. However, the 
magnitude of the differential pressure cell output, AP in 
equation (51) , dominates the value of computed level so the 
resulting effect on generator level of a reactor building 
temperature change of this magnitude would be only .7 
percent. A steam leak into the reactor building would 
significantly raise the ambient humidity, making all 
instrumentation suspect of erratic operation and place the 
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plant at risk of uncovering the core . Under the 
circumstances, a .7 percent change in indicated steam 
generator level would not be as important as the other 
effects on the reactor plant. 

Obviously, none of the simplifications and assumptions made 
by the conventional instrument envisioned for comparison 
with the smart instrument produce an error in excess of 1.6 
percent. This small deviation from the "true" steam 
generator level is really too small for the operators to 
notice, so the smart instrument measuring steam generator 
level does not provide a realistic increase in indication 
accuracy. It still would be useful for warning the operator 
of the approach to a flooding situation in the drain pipe. 

6 . 2 Smart Instrument Design 

This research indicates there are several considerations to 
observe when designing a smart instrument: 

•select the variables for measurement; 

•choose a model that incorporates as much detail about 
the physical process as necessary to achieve the 
measurement accuracy desired and the computational 
resources of available microprocessors; 

•evaluate the model stability and accuracy of results; 

•compare the smart instrument with a conventional 
instrument to evaluate the magnitude of the accuracy 
improvement ; 

•compare information available to the operator using 
the conventional instrument to that provided by the 
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smart instrument, even if not significantly more 
accurate, the designer may still opt for the smart 
instrument because of the extra information it 

6 . 3 Flooding 

The suggested model describing steam generator level 
dynamics, to the extent it closely approximates the actual 
variables involved, is capable of calculating the liquid 
level in the drain pipe between the steam generator and the 
condensing pot directly. By observing when the liquid level 
reaches or almost reaches the top of the drain pipe, the 
smart level instrument envisioned by this research could 
warn the operator that its level output data is on the verge 
of unreliability. It is important to note, that the ability 
of the smart instrument to compute a solution of equation 
(49) for the state variables is predicated on the condition 
that flooding in the drain pipe does not occur, since 
equations (28) and (29) , used to relate the liquid height in 
the drain pipe and the pressure gradient between the 
condensing pot and steam generator to the superficial 
velocities in the pipe, are not valid after flooding occurs. 
The pertinent question then is how does the smart instrument 
recover from just such a flooding condition since its model 
is now no longer valid? 

To recover from a flooding condition, the smart instrument 
should be programmed with a routine to estimate the steam 
generator water level based on the reduced reference leg 
height discussed in section 4.1 and depicted in figure 5. 

The instrument should warn the operator that flooding has 
occurred, but still be able to present its best estimate for 
water level. Next, the routine used to find the superficial 
velocities from the tables that have been pre-computed must 
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be capable of responding to the case when the input 
arguments, liquid height in the drain pipe and pressure 
gradient, fall outside the bounds of the tabular data. The 
smart instrument should then try to calculate level after 
discarding its prior data, in effect re-zeroing itself. It 
should continue to zero itself and try to calculate 
transient values of the state variables until the a the 
drain pipe liquid level is reduced and the pipe is no longer 
flooded. 

Conclusions 



This research discusses how and why smart instruments might 
be applied to nuclear power plants. The features that 
distinguish a smart instrument from a conventional 
instrument are a process model and the ability of the smart 
instrument to make decisions based on plant variables and 
model output. A good model of the physical process enables 
the instrument to monitor measurement results and warn the 
user if the data received is not consistent with redundant 
sensors or trends indicated by sensors for different 
variables. The important issues to consider when designing 
a smart instrument were explored and actually illustrated by 
a specific example: design of a smart instrument for steam 

generator water level measurement. A model of steam 
generator level variations as a function of power was 
proposed to reduce the number of assumptions necessary for 
level computation and to infer unmeasured variables germane 
to the reliability of the remaining assumptions. The Prism 
program was used to generate simulated level data, from 
which an estimated differential pressure cell output signal 
was extracted. The differential pressure cell signal and 
other Prism data were input to the smart level instrument. 
The "smart" level produced was compared to a conventional 
level instrument's output. The smart level signal was only 
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marginally more accurate than the conventional level signal 
but the smart instrument is capable of inferring important 
variables that the conventional instrument cannot. Finally, 
the smart instrument must be capable of recovering from 
brief operation outside the limits of its model, otherwise 
the designer risks that certain transients could render the 
instrument's subsequent output unreliable. 
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Chapter 7 



Recommendations for Future Work 



There has not been sufficient time during the term of this 
research to pursue all the issues involved in smart 
instrumentation implementation as thoroughly as desired. 
Topics related to a smart level instrument worthy of 
additional investigation include: determining the criteria 

affecting the theoretical stability of the one-step, 
implicit difference equation (equation (48)) used to compute 
the solution to the linearized equations, investigating 
model behavior at lower steam generator pressures, 
numerically analyzing of the solution method, and obtaining 
experimental data on flow transitions and pipe liquid height 
for two phase counter-current flow and comparing the 
experimental results to the predicted results. Applying the 
design methodology suggested in this research or another 
method to the measurement and display of other nuclear 
reactor plant variables, such as reactivity measurement and 
control, would be of interest for the continued study of 
smart instruments in general. 

7 . 1 Theoretical Stability 

Before implementing the proposed model for steam generator 
level measurement, the specific conditions for stability of 
equation (24) , the general differential equation, and 
equation (48) the one-step implicit difference, should be 
determined . 

7 . 2 Operation at Lower Pressures 
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The steam generator level model developed in this research 
was only applied to level measurement at normal steam 
generator operating pressure, about 7.0 to 7.6 MPa. For the 
smart instrument to really replace conventional level 
instruments, it must be capable of operating over the entire 
range of generator pressures encountered in practice, 
extending down to atmospheric pressure. To investigate 
model performance at other pressures, solutions to the 
momentum balance equations, (28) and (29) , must be computed 
for several pressures between 7 . 6 MPa and atmospheric since 
the liquid and vapor densities, taken as constants when 
iteratively solving equations (28) and (29) , are strong 
functions of pressure in the steam generator. 

7 . 3 Physical Testing 

As the comparison of results in Chapter 5 illustrated, the 
smart level instrument employing a detailed model was only 
marginally more accurate than a conventional level 
instrument using a simplified model. Still, the smart 
instrument is capable of calculating the transient liquid 
level in the drain pipe and can therefore detect the onset 
of flooding in the drain pipe and subsequent erroneous level 
indication. Since the ability to predict the drain pipe 
liquid level is the distinguishing feature of the smart 
instrument, the accuracy of the prediction should be 
evaluated by physical experiment. The physical data 
obtained could be used to improve the model used as well. 

7 . 4 Numerical Analysis 

Computing the solution to equation (48) is very time 
consuming. The solution for the positive pressure ramps 
requires several hours to complete on a DECstation 3100, a 
particularly fast workstation. It would be potentially 
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beneficial to further analyze the equations involved and 
construct an algorithm designed to enhance the speed of the 
solution. Certainly, the speed of the calculations could be 
significantly enhanced by the use of a compiled high level 
computer language for the computations. 

7 . 5 Application to Other Variables 

The application of smart instruments to measurement and 
control of other nuclear power plant variables should be 
investigated. In addition to evaluating and displaying the 
value of a simple plant variable, use with complicated 
systems like reactivity or steam generator water level 
control should be explored. A smart instrument could be a 
key component of an expert system used for power plant 
operation and casualty control. The smart instrument could 
be responsible for redundancy and consistency checks while 
an external computer evaluates the over all system, 
postulates probable outcomes to certain actions, and can 
recommend corrective action. Specific applications 
utilizing the ability of smart instruments to communicate 
with each other should also be explored. Sensors could be 
networked together and either polled by a control system or 
signal the control system when they have data to transmit. 
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Appendix 



The PRO-MATLAB program, numerical processing software based 
on LINPACK and EISPACK codes, is used to solve equation (48) 
for each of the postulated transients. PRO-MATLAB is an 
optimized program that is well-suited for high speed 
calculations on large matrices and vectors . The program is 
capable of processing commands that have been saved to a 
text file, called a macro, and allows the user to create his 
own functions. Neither functions or macros are compiled, 
but they were used to develop prototype algorithms to 
compute solutions to equation (48) . The same routines could 
be compiled and used by the smart level instrument to 
perform similar calculations. The matlab macros and 
functions used are included in the following pages. 
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%macro to start calculating a solution to equation (48) 
stmdat 

%gets rhol_list , rhov_list , Hl_list, Hv_list, pressure 
variables 

%gets table_h and table_j?ress values (for jv and jl data) , 
jvdat and jldat 
dt=l* . 5; 

L=2 . 5; %length of drain pipe 

Dt= . 0254 ; %drain pipe dia 

Dc=2* . 0254 ; %condensing pot dia 

Lc=3 . 5*2 . 54e-2; %length of condensing pot 

dpdx=-73 .27493046894671; 

startpsg=7 . 34e6; 

endpsg=7 . 0400e6 ; % 

stoptrans=225/dt ; 

minutes=4 ; %number of minutes to simulate 
input= [startpsg+dpdx*L; . 71227377; 0 . 00546943] ; 

deltas= [-30; . 001; . 0005] ; %deltas to use in computing 

derivatives 

X_= [ ] ; 

f _= [ ] ; 
prsg= [ ] ; 

%X is the set of intermediate values of the result for each 
time %step 

for n=-5/dt : minutes* 60/dt 
%controls psg(t)-> 
if n<=0 

psg=startpsg; %initial value 
elseif n<stoptrans+l 

psg=startpsg-n* (startpsg-endpsg) /stoptrans; %for ea time 
%step during power chg 
else 

psg=endpsg; 

end 

X_= [X_ input ] ; 
prsg= [prsg; psg]; 
psg 

% display dp/dx and dp/dx + dp 
(input (1) -psg) /L 
(input (1) -psg) /L+deltas (1) 

[X, f ] =solveit (psg, 70000, dt, input, L, Dc, Dt, Lc, deltas, . . . 
jvdat, jldat, rhol_list, rhov_list , Hl_list , Hv_list , . . . 
table_press, table_h, pressure) ; 
input =X; 

disp (' finished time step') 
n*dt 

f_=[f_ f]; 
end 
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%this is the macro stmdat 

%get steam table data in flat file form 

load stmdata.prn 

n=16; %number of pressure data pairs 
pressure=stmdata ( 1 : n) ; 

rhol_list=ones (1 : n) ' . /stmdata (n+1 : 2*n) ; 
rhov_list=ones (1 :n ) ' . /stmdata (2*n+l : 3*n) ; 

Hl_list=stmdata (3*n+l : 4*n) ; 

Hv_list=stmdata (4*n+l : 5*n) / 

%load outdata3.prn 
load outdata4.prn 

nl=41 ; %number of pressure data pairs 

table_h= [ 0 ; .001; .0015; .002; .0025; .003; .0035; .004; .008; . . 
. 012; . 016; . 020; . 023; 025] ; %order of liquid heights 

table_press= [600; 570; 540; 510; 480; 450; 420; 390; 360; 330; 300 
270; 240; 210; 180; 150; 120; 90; 60 ; 30; 0; -30; -60 ; -90 . 1;-120 .1; 
-150 . 1 ; -180 . 2 ; -210 . 2 ; -240 . 2; -270 . 2; -300 . 3; -330 . 3 ; -360 . 3 ; 
-390. 3; -420. 4; -450. 4; -480. 4; -510. 4; -540. 5; -570. 5; -600. 5] 

jvdat=outdata4 (1 : 14*nl) ; 
jldat=outdata4 (14*nl+l : 28*nl) ; 
end 
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function [X, f ] =solveit (psg, dpsg, dt , inputx, L, Dc, Dt , Lc, . . . 
d, jvdat, jldat,rhol_list,rhov_list,Hl_list,Hv_list, . . . 
tablejpress, table_h, pressure) ; 

%function to solve the equation and compile the answers 

%psg is initial generator pressure 
%dpsg is delta generator pressure 
%dt is time interval of manuever 
% inputx is inputx= [pc; alpha_c;h] 

%L is length of drain pipe 
%Dc is condensing pot diameter 
%Dt is drain pipe diameter 

%d is deltas of x variables for numerical differentiation 
%d= [ dpc ; dalpha_c ; dh ] [-30; .0005;. 0005] 

%you get jvdat, jldat, rhol_list, rhov_list, Hl_list, 
Hv_list , table_press , 

%table_h from running stmdat 

%pressure is the list of pressures for which steam table 
data is available in 

%the file stmdat a. prn read by stmdat. m 

%define given deltas to use for numerical derivatives 
dp=d ( 1 ) ; 
da=d (2) ; 
dh=d (3) ; 

% 

%for n=l : 5 

% psg=psg-n*dpsg/ 600 ; %f or each time step 
%define state variables 
alpha_c=inputx (2) ; 
pc=inputx (1) ; 
hl=inputx (3) ; 

%now build the df/dx vector 

df ll=f irst (alpha_c, psg, pc, hi, L, Dc, Dt , dp, jvdat, jldat , . . . 
table_press, table_h) ; 

df 12=second (alpha_c, psg, pc, hi, L, Dc, Dt, da, jvdat , jldat, . . . 
table_j?ress, table_h) ; 

dfl3=third(alpha_c,psg,pc,hl,L,Dc,Dt,dh, jvdat, jldat, . . . 
table_j?ress, table_h) ; 

df21=fourth (alpha_c, psg, pc, hi, L, Dc, Dt, dp, jvdat, jldat, . . . 
table_press, table_h) ; 

df 22=f ifth (alpha_c, psg, pc, hi, L, Dc, Dt , da, jvdat, jldat , . . . 
table_press, table_h) ; 

df 23=sixth (alpha_c, psg, pc, hi, L, Dc, Dt , dh, jvdat , jldat , . . . 
table_press, table_h) ; 

df31=seventh(alpha_c,psg,pc,hl,L,Dc,Dt,dp, jvdat, jldat, . . . 
table_press, table_h) ; 

df 32=eighth (alpha_c, psg, pc, hi, L, Dc, Dt , da, jvdat , jldat , . . . 
table^press, table_h) ; 

df 33=ninth (alpha_c, psg, pc, hi, L, Dc, Dt , dh, jvdat , jldat , . . . 
table_press, table_h) ; 
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df= [df 11 df 12 df 13/ df21 df22 df23;df31 df32 df33]; 

% 

%get the f vector 

f=f indf (alpha_c, psg, pc, hi, Dc, Dt , L, jvdat , jldat , . . . 
table^press, table_h) 

% 

%get the A vector 
% 

A=buildA (pc, alpha_c, hi, rhol_list, rhov_list , Hl_list, . 
Hv_list, Dc, Dt, L, Lc, pressure) ; 

X=inv((A/dt - df) ) *f + inputx; 

disp (' finished time step') 

end 
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function [df ldpc] =f irst 

(alpha_c, psg, pc, hi, L, Dc, Dt , dp, j vdat , jldat , press , givenh) 

% 

%function to compute the derivative of fl wrt pc 

% 

%f unctions first.m through ninth. m use similar arguments: 

% 

%alpha_c is the vapor volume fraction in the condensing pot 

%psg is the steam generator pressure 

%pc is the condensing pot pressure 

%hl is the liquid height in the drain pipe 

%L is the length of the drain pipe 

%Dc is the diameter of the condensing pot (a horiz cylinder) 

%Dt is the diameter of the drain pipe 

%dp or da or dh are the deltas used to compute the 

derivatives %wrt 

%pc,alpha_c, or hi 

%press and givenh are the P' and h data used to compute the 

jvdat %and jldat 

%tables 



theta= . 0164 ; 

% 

%get steam data 

[a, b, Tsat , mu, k, rhovl, rhol] = stminfo (pc); 

[a,b, Tsat,mu, k, rhov2, rholj = stminfo (pc+dp) ; 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt , pc, theta) ; 
mdotw2 = mdotwl; %same as 

mdotw2=f indmdot (alpha_c, Dc, Dt, pc+dp, theta) ; 

% 

%find jv 
% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

[jl2 jv2 ] = 

findjvjl (hi, (pc-psg) /L+dp, givenh, press, jvdat, jldat) ; 

% 

A= (Dt/2) A 2*pi; %area of pipe 
%f l=rhovl * jvl *A-mdotwl 

dfldpc = ( (rhov2* jv2*A - mdotw2) - (rhovl* jvl *A - mdotwl)) /dp; 
end 
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function [dfldac] =second (alpha_c, psg, pc, hi , L, Dc, . . . 

Dt, da, jvdat , jldat, press, givenh) 

% 

%function to compute the derivative of fl wrt alpha_c 
% 

%function [dfldac] = second 

(alpha_c , psg, pc, hi, Dc, Dt , da, jvdat, jldat ) 

% 

theta= . 0164 ; 



%get steam data 

[a, b, Tsat , mu, k, rhovl, rhol] = stminfo (pc); 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 
mdotw2 = findmdot (alpha__c+da, Dc, Dt, pc, theta) ; 

% 

%find jv 
% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press , jvdat, jldat) ; 

[ j 12 jv2 ] = 

findjvjl (hi, (pc-psg) /L, givenh, press , jvdat, jldat) ; 

% 

A= (Dt/2) A 2*pi; %area of pipe 
%f l=rhovl* jvl*A-mdotwl ; 

dfldac = ( (rhovl* jv2*A - mdotw2) - (rhovl* jvl *A - mdotwl) ) /da; 
end 
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function [dfldh] =third (alpha_c, psg, pc, hi , L, Dc, . . . 

Dt, dh, jvdat , jldat, press, givenh) 

% 

%function to compute the derivative of fl wrt h 

% 

%function [dfldpc] = third 

(alpha_c, psg, pc, hi, L, Dc, Dt , dh, jvdat, jldat) 

theta= . 0164 ; 

% 

%get steam data 

[a, b, Tsat , mu, k, rhovl , rhol] = stminfo (pc); 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 

% 

%find jv 
% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

[ jl2 jv2] = 

findjvjl (hl+dh, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

% 

A= (Dt/2) /s 2*pi; %area of pipe 

dfldh = ( (rhovl* jv2*A - mdotwl) - (rhovl* jvl *A - mdotwl) ) /dh; 
end 
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function [df2dpc] =fourth (alpha_c, psg, pc, hi, L, Dc, Dt, . . . 
dp, jvdat, jldat , press, givenh) 

% 

%function to compute the derivative of f2 wrt pc 
% 

%function [dfldpc] = fourth 

(alpha_c, psg, pc, hi, Dc, Dt , dp, jvdat, jldat) 

theta=. 0164; 

% 

%get steam data 

[Hll , Hvl , Tsat , mu, k, rhovl , rhol] = stminfo (pc); 

[H12 , Hv2 , Tsat , mu, k, rhov2 , rhol] = stminfo (pc+dp) ; 

% 

%find mdot 

% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 
mdotw2 = mdotwl ;% same as 

mdotw2=f indmdot (alpha_c, Dc, Dt, pc+dp, theta) ; 

% 

%find jv 
% 

[jll jvl] = 

f ind jvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

[ j 12 jv2] = 

f indjvjl (hi, (pc-psg) /L+dp, givenh, press, jvdat, jldat) ; 

% 

%f indQ 
% 

Ql=f indQ (pc) ; 

Q2=findQ (pc+dp) ; 

A= (Dt/2) /v 2*pi; %area of pipe 
f 2=rhovl* jvl*Hvl*A-mdotwl*Hll— Q1 ; 
df2dpc = 

( (rhov2* jv2*Hv2*A-mdotw2*Hl2-Q2) - (rhovl* jvl*Hvl*A-mdotwl*Hll 

Ql) ) /dp; 
end 
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function [df2dac] =fifth (alpha_c, psg, pc, hi, L, Dc, Dt, . . . 
da, jvdat , jldat, press, givenh) 

% 

%function to compute the derivative of f2 wrt alpha_c 

% 



theta=. 0164 ; 

% 

%get steam data 

[Hll , Hvl , Tsat , mu, k, rhovl , rhol] = stminfo (pc); 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 
mdotw2 = findmdot (alpha_c+da, Dc, Dt, pc, theta) ; 

% 

%find jv 

% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

% 

%f indQ 

% 

Ql=f indQ (pc) ; 

A= (Dt/2) A 2*pi; %area of pipe 
% f 2=rhovl * jvl *Hvl *A-mdotwl *H1 1 -Q1 
df2dac = 

( (rhovl* jvl *Hvl*A-mdotw2*Hll-Ql) - (rhovl* jvl*Hvl*A-mdotwl*Hll 

Ql) ) /da; 
end 
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function [df2dh] =sixth (alpha_c, psg, pc, hi, L, Dc, Dt, . . . 
dh, jvdat, jldat, press, givenh) 

% 

%function to compute the derivative of f2 wrt h 
% 

%function [dfldh] = third 
(alpha_c,psg, pc, hi, Dc, Dt, dh, jvdat, jldat) 

theta= . 0164 ; 

% 

%get steam data 

[Hll , Hvl , Tsat , mu, k, rhovl , rhol] = stminfo (pc); 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 

% 

%find jv 

% 

[jll jvl] = 

f indjvjl (hi, (pc-psg) /L, givenh, press , jvdat, jldat) ; 

[ j 12 jv2] = 

findjvjl (hl+dh, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

% 

%f indQ 

% 

Ql=f indQ (pc) ; 

A= (Dt/2) / '2*pi; %area of pipe 
f 2=rhovl* jvl*Hvl*A-mdotwl*Hll-Ql; 
df 2dh = 

( (rhovl* jv2*Hvl*A-mdotwl*Hll-Ql ) - (rhovl* jvl*Hvl*A-mdotwl*Hll 

-Ql) ) /dh; 

end 
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function [df3dpc] =seventh (alpha_c, psg, pc, hi, L, Dc, Dt, . . . 
dp, jvdat , jldat, press, givenh) 

% 

%function to compute the derivative of f3 wrt pc 
% 



theta=. 0164; 

% 



%get steam data 

[a, b, Tsat,mu, k, rhovl, rholl] = stminfo (pc); 

[a, b, Tsat , mu, k, rhov2, rhol2 ] = stminfo (pc+dp) ; 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 
mdotw2 = findmdot (alpha_c, Dc, Dt, pc+dp, theta) ; 

% 

%find jv 

% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press , jvdat, jldat) ; 

[jl2 jv2 ] = 

findjvjl (hi, (pc-psg) /L+dp, givenh, press, jvdat, jldat) ; 

% 

A= (Dt/2) A 2*pi; %area of pipe 
f 3=rhol 1 * j 11 *A+mdotwl ; 

df3dpc = ( (rhol2* jl2*A + mdotw2) - (rholl* jll*A + mdotwl) ) /dp; 
end 
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function [df 3dac] =eighth (alpha_c, psg, pc, hi, L, Dc, Dt, da, . . . 
jvdat, jldat, press, givenh) 

% 

%function to compute the derivative of f3 wrt alpha_c 
% 

theta= .0164; 

% 

%get steam data 

[a,b,Tsat,mu,k,rhovl, rholl] = stminfo (pc); 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 
mdotw2 = findmdot (alpha_c+da, Dc, Dt, pc, theta) ; 

% 

%find jv 
% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

[ j 12 jv2] = 

findjvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 

% 

A= (Dt/2) / '2*pi; %area of pipe 
%f 3=rholl* jll*A+mdotwl ; 

df3dac = ( (rholl* jl2*A + mdotw2) - (rholl* jll*A + mdotwl) ) /da; 
end 
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function [df3dh] =ninth (alpha_c, psg, pc, hi, L, Dc, Dt, . . . 
dh, jvdat, jldat, press, givenh) 

% 

%function to compute the derivative of f3 wrt h 
% 

theta=. 0164; 

% 

%get steam data 

[a, b, Tsat, mu, k, rhovl , rholl] = stminfo (pc); 

% 

%find mdot 
% 

mdotwl = findmdot (alpha_c, Dc, Dt, pc, theta) ; 

% 

%find jv 
% 

[jll jvl] = 

findjvjl (hi, (pc-psg) /L, givenh, press, jvdat , jldat) ; 

[ jl2 jv2 ] = 

findjvjl (hl+dh, (pc-psg) / L, givenh, press, jvdat, jldat) ; 

% 

A= (Dt/2) A 2*pi; %area of pipe 

df 3dh = ( (rholl* jl2*A + mdotwl) - (rholl* jll*A + mdotwl) ) /dh; 
end 
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function [f] =findf (alpha_c, psg, pc, hi , Dc, Dt , . . . 

L, jvdat , jldat , press , givenh) 

%function to compute the f vector 
% 

%alpha_c is the vapor volume fraction in the condensing pot 
%psg is steam generator pressure 
%pc is condensing pot pressure 
%hl is liquid height in drain pipe 

%Dc is condensing pot diameter (a horiz cylinder) 

%Dt is drain line diameter 
%L is length of drain line 

%press and givenh are the P' and h data used to compute the 

jvdat %and jldat 

%tables 



^constants 
theta= .0146; 

% 

%get steam data 

[HI, Hv, Tsat , mu, k, rho_v, rho_l] =stminfo (pc) ; 

% 

%find mdot 

% 

mdotw = findmdot (alpha_c, Dc, Dt, pc, theta) ; 

% 

%find jv, jl 

[jl j v l = findjvjl (hi, (pc-psg) /L, givenh, press, jvdat, jldat) ; 
% 

%find Q 
Q=findQ (pc) ; 

% 

A= (Dt/2) A 2*pi; %area of pipe 
% 

%f ind fl 

f l=rho_v* jv*A-mdotw; 

% 

%find f2 

f 2=rho_v* jv*Hv*A-mdotw*Hl-Q; 

% 

%find f3 

f 3=rho_l * j 1 *A+mdotw; 

% 

f=[fl;f2;f3] ; 
end 
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function [A] = buildA (pc, alpha_c, hi , rhol_list , rhov_list , . . . 
Hl_list , Hv_list, Dc, Dt , Lt , Lc, pssr) 

% 

%function that builds the A matrix, xxx list are the stm 
table %lists of 

%data for the pressures read in from stmdata.prn by stmdat.m 
%pssr is the range of pressures used to extract steam table 
data 

% 

%pc is pressure in condensing pot 

%alpha_c is vapor volume fraction in condensing pot 
%hl is liquid hight in the drain pipe 
%Lt is the length of the drain pipe 

%Lc is the length of the condensing pot (assumed to be a 
horiz %cylinder) 

% 

%uses : 1. stminfo.m 

% 2. deriv.m 

% 

Vc= (Dc/2) ^2*pi*Lc; %volume of condensing pot 
[Hi Hv Tsat mu k rhov rhol] =stminfo (pc) ; %stm table info 

% 

%note : deriv (pssr, rhol_list,pc, . 01) computes the partial 
deriv of %rhol wrt pc at pc 

% compute all 

all= ( (l-alpha_c) *deriv (pssr, rhol_list , pc, .01) +alpha_c* . . . 
(deriv (pssr, rhov_list , pc, . 01) ) ) *Vc; 

% 

%compute al2 
al2= (rhov-rhol) *Vc; 

% 

% compute a2l 

a21= ( (l-alpha_c) *deriv (pssr, rhol_list . *Hl_list , pc, . 01 ) + . . . 
alpha_c* (deriv (pssr, rhov_list . *Hv_list,pc, . 01) ) -1) *Vc; 

% 

% compute a22 

a22= (rhov*Hv - rhol*Hl) *Vc; 

% 

%compute a33 

a33=Dt/2*Lt*rhol* (sin(acos (l-2*hl/Dt) ) ) "2/sqrt (hl/Dt* (1-hl/D 
t) ) ; 

% 

%complete A 

A= [all, al2, 0 ; a21, a22, 0; 0, 0, a33] ; 
end 
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function y = deriv (x, y, input , delta) 

% 

^computes the derivative of y wrt x at the input value of x, 
%using an x step 
%of delta 

%format is y = deriv (x, y, input , delta) 

y = (spline (x, y, input+delta) -spline (x, y, input) ) /delta; 

end 
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function [hi, hv, Tsat , mu, k, rho_v, rho 1] = stminfo (p) 

% 

%function to compute the properties of water given pressure 

% 

%for pressure in the range 1.5MPa to 9Mpa 
% 

%function [hi, hv, Tsat , mu, k, rho_v, rho_l] = stminfo (p) 

k=62 6e— 6 ; %thermal conductivity 
%steam table information 

press = [1 . 5e6; 2e6; 2 . 5e6; 3e6; 3 . 5e6; 4e6; 4 . 5e6; 5e6; 5 . 5e6; . . . 
6e6; 6 . 5e6; 7e6; 7 . 5e6; 8e6; 8 .5e6; 9e6] ; 

rhol= [866. 6262; 849. 9065 ; 835. 4219; 822. 1656; 810. 044 6; . . . 

7 98. 658 3 ; 787. 9 600 ; 777. 7259 ,*7 67.8722 ,-758.322 6; 74 9. 0 637; . . . 

74 0.0281 ;731.1545 ; 722. 4390 ; 7 13. 7759 ,-705.2684] ; 

rhov = [7 .5953; 10 . 0462; 12.6855; 15.0083; 17.5346;... 
20.1005; 22.7066; 25.3614; 28.0662; 30.8261; 33.6485;... 
36.5323; 39.4836; 42.5080; 45.4339; 48.7924]; 

Hl=1.0e+03*[. 8447; 0.908 6; 0.9618; 1.0 08 4; 1.04 98; 1.0874; . . . 
1.1221;1.1545;1. 1849; 1.2137; 1.2417; 1.2674; 1.2 927; 1.3171; . . . 

1 . 3404; 1 . 3637] ; 

Hv=1.0e+03* [2. 78 99; 2. 7 972; 2. 800 9; 2. 8023; 2. 8 020; 2. 8003; . . . 

2. 7977; 2. 7942; 2. 7899; 2. 7850; 2. 7795; 2. 7735; 2. 7669; 2. 7599; . . . 
2. 7523 ,*2.7446] ; 

TSAT = [198. 29; 212. 37; 223. 91; 233. 84; 242. 54 ;250.33; 257. 41; . . . 
263. 91; 2 69. 93; 275. 55; 280. 82; 285. 7 9; 290. 5 ,*294.97 ,*299.14; 303. 3 
1]; 

%interpolate to get precise T and P in generator for given 
power % level 

%i and i+1 are the indices in the list that bracket the 
value % sought 
i=comp (p, press) ; 
if i==length (press) 
rho_v=rhov (i) ; 
rho_l=rhol (i) ; 
hv=Hv ( i ) ; 
hl=Hl (i) ; 

Tsat=TSAT (i) ; 
else 

rho_v= (p-press (i) ) . / (press (i+1) -press (i) ) . * (rhov (i+1) — . . . 
rhov (i) ) +rhov (i) ; 

rho_l= (p-press (i) ) . / (press (i+1) -press (i) ) . * (rhol (i+1) — . . . 
rhol (i) ) +rhol (i) ; 

hv= (p-press (i) ) . / (press (i+1) -press (i) ) . * (Hv (i+1) -Hv (i) ) + . . . 
Hv (i) ; 

hl= (p-press (i) ) . / (press (i+1) -press (i) ) . * (HI (i+1) -HI (i) ) +. . . 
HI (i) ; 
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Tsat= (p-press (i) ) . / (press (i+1) -press (i) ) . * (TSAT (i+1) - . . . 
TSAT (i) ) +TSAT (i) ; 
end 

mu=visc (Tsat, rho_l) ; 
end 
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function [index] =comp (value, list) 

%function index=compare (value, list) 

%function that compares each element of value with all the 
values in a 

%sorted list and returns the a vector containing the index 
of the closest 

^number below each element of value in list 

index=zeros ( 1 , length (value) ) ; 
for j=l : length (value) ; 

for i=l : length (list ) ; 
if value (j) >=list (i) 
index ( j ) =i ; 
end 
end 

if index (j)==0 
disp (' ' ) 

disp(' number is not in given range') 
end 

end 

end 
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function [Q]=findQ(p) 

%function to compute heat transfer rate as a function of 
pssr in %condensing pot 
%function [Q]=findQ(p) 

%constants 
g=9 . 807 ; 

Tamb=50 ; 

%wtd length. length is 3.5 inches half of which is 1.75 
inches=5 . 08e-2 
1=3. 5/2*2. 54e-2; 

%l=3.81e-2; 

%rin=7 . 42e-2 ; 

%rout=7 . 62e-2; 

rin=l . 939*2 . 54e-2/2 ; %based on sched 80, 2 in diameter 
rout=2 . 157*2 . 54e-2/2; 
ksteel=50; 
hout=ll ; 

% 

%get steam data 

[hi, hv, Tsat,mu, k, rho_v, rho_l] =stminfo (p) ; 

% 

%hfg=hv-hl; 

%compute hbar 

%Nu= (hfg . *rho_l . A 2*g*l / '3 . / (4*mu*k. * (Tsat-Tsurf ) ) ) . ^ . 25; 
%hbar=4/3*k*le3/l*Nu; %le3 needed since k is in kw 
% compute UA 

%UA= ( (1/ (2*pi*rin*l*hbar) ) + (log (rout/ rin) / (2*pi*l*ksteel) ) + . 

(1/ (2*pi*rout*l* . . . 

%hout ) ) ) * (-1) ; 

UA= ( (log (rout/rin) / (2*pi*l*ksteel) ) + (1/ (2*pi*rout*l* . . . 

hout)))*(-l); 

% compute Q 

Q=UA. * (Tsat-Tamb) ; %mass flow rate 
end 
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function [mdotw] = findmdot (alpha_c, Dc, Dt, pc, theta) 

% 

%function to compute the mass flow rate in the weir 
% 

%alpha_c is the vapor volume fraction in the condensing pot 

%Dc is the diameter of the condensing pot 

%Dt is the diameter of the drain pipe 

%pc is the pressure in the condensing pot 

%theta is the elevation angle of the drain pipe 
% 

%uses : 

% l.stminfo.m 

% 2.interp.m 

%constants 
g= 9.807; 
a 1=2 ; 

%find the liquid density 

[a, b, c, d, e, f , rhol] = stminfo (pc) ; 

%find the phi_c that corresponds to the given alpha_c 
phi_c=interp (alpha_c) / %call special routine to find phi_c 



hlc=Dc/2* (1-cos (phi_c) ); %liq ht in condensate pot 
% 

if hlc> (Dc-Dt) /2 & hlc< (Dc+Dt) /2 %make sure hlc in correct 
range 

hw=hlc- (Dc-Dt) / 2; 
else 

hw=-l;%tube is plugged or empty 
end 
% 

if hw~=-l %tube not plugged 
phi_w = acos (l-2*hw/Dt) / 

alpha_b = l-(phi_w - sin (phi_w) *cos (phi_w) ) /pi ; 
else 

alpha_b = 1; 
end 

%alpha_b is the vapor volume fraction at the weir 

% 

%use eqn () 
mdotw = 

rhol* (Dt/2) A 2*pi* (l-alpha_b) *sqrt (g*Dt*cos (theta) /al) / 
end 
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function [phi_c] =interp (alpha_c) 

% 

%function [phi_c] =interp (alpha_c) 

% 

%function that finds phi given alpha using bisection search 

phi_c=pi/2 ; 
bottom=0; 
top=pi ; 
dif f=l ; 

while abs (diff ) >5e-9, 

diff=alpha_c- (1- (phi_c-sin (phi_c) *cos (phi_c) ) /pi) ; 
if diff < 0 
bott om=phi_c ; 

phi_c=phi_c+ (top-phi_c) /2; 
else 

top=phi_c; 

phi_c=phi_c- (phi_c-bottom) /2; 
end 
end 
end 
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function [jl, jv]— fiadjvjl (hi, Pprime, h, press, jvdat, jldat) 

% 

%This macro does a 2 way interpolation on jv and jl data 
compiled %as a function of both hi and P' 

%hl is liquid height in the drain pipe 
%Pprime is dP/dx for the tube 

%press and givenh are the P' and h data used to compute the 
jvdat %and jldat tables 
% 

%function [jl,jv]= findjvjl (hi, Pprime, h, press, jvdat , jldat) 

% 

%find out where hi and Pprime fall in the range of h' s and 
%press's available 
% 

%note that h and press contain all the possible values of hi 
and %Pprime 
% 

a=find(hl<=h & hl>=h(l)); %first element of a is the first h 
%higher than given hi 

b=f ind (Pprime>=press & Pprime<=press (1) ) ; 

% 

%now select the 4 jv/jl values that bracket the hi and P' 
given 
% 

if length (a) ~=0 & length (b) ~= 0 %hl and Pprime were in the 
%correct ranges 
% 

%fix pressure index (at lower bracket value, one less than 
b (1) 

if b (1) ==1 

b_=l; %in case pressure index already at lowest 

% value, don't decrement 
else 

b_=b (1) — 1; %fix pressure index at next lower value 
end 

%drop hi index to lower bracket 
if a (1) ==1 

a_=a(l)+l; %so a_-2 will be 0,not negative 
else 

a_=a ( 1 ) ; 
end 

jvhl=jvdat (41* (a_-2) +b_) ; %lower bracket hi and Pprime data 
jlhl= jldat (41* (a_-2) +b_) ; 

%bump hi index back to upper bracket 
if b (1) ==1 

b_=l; %in case pressure index already at lowest 

%value, don't decrement 
else 

b_=b (1) -1; %keep pressure index at next lower value 
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end 

a_=a (1); 

jvh2=jvdat (41* (a_-l) +b_) ; %upper bracket hi, lower bracket 
%Pprime data 

jlh2= jldat (41* (a_-l) +b_) / 

%bump pressure to upper bracket 
% 

%fix pressure index (at upper bracket) 
if b (1) ==1 

f_=2; %in case pressure index already at lowest 

% value, bump up 
else 

f_=b(l); %put pressure index at upper bracket 
end 

%drop hi index to lower bracket 
if a(l)==l 
g_=a (1) +1/ 
else 

g_=a (1) ; 
end 

jvh3= jvdat (41* (g_-2) +f_) ; %lower bracket hi, upper bracket 
P' 

jlh3= jldat (41* (g_-2) +f_) ; 

%bump hi index back to upper bracket 
if b (1) ==1 

f_=2; %in case pressure index already at lowest 

% value, bump up 
else 

f_=b(l); %put pressure index at upper bracket 
end 

g_=a (1); 

jvh4= jvdat (41* (g_-l ) +f_) ; %upper bracket hi and P' 
jlh4= jldat (41* (g_-l) +f_) ; 

% jvh4= jvdat (11* (g_—l) +f_) ; %upper bracket hi and P' 

% jlh4=jldat (11* (g_-l) +f_) ; 

%interpolate between h values 
if a ( 1 ) ==1 
hl=h (1) ; 
else 

hl=h (a (1) -1) ; 
end 

h2=h (a (1) ) ; 
h3=hl; 
h4=h2 ; 

%these if' s are required since hi might = h2 and the same 
for h3 %and h4 
if h2==hl 
jv5= jvhl; 
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else 

jv5= (hl-hl) . / (h2-hl) * ( jvh2-jvhl) + jvhl; 
end 

if h2==hl 
jv6= jvh3 ; 
else 

jv6=(hl-h3) . / (h4-h3) * ( jvh4-jvh3) +jvh3; 
end 

if h2==hl 
jl5= jlhl ; 
else 

jl5= (hl-hl) . / (h2-hl) * ( jlh2-jlhl) + jlhl; 
end 

if h2==hl 
jl6= jlh3 ; 
else 

jl6= (hl-h3) . / (h4-h3) * ( jlh4-jlh3) + jlh3; 
end 

%interpolate between 5 6 for Pprime 
if b (1) ==1 

pl=press (b (1) ) ; 
else 

pl=press (b (1) -1) ; 
end 

if b (1) ==1 

p2=press (b (1) +1) ; 
else 

p2=press (b (1) ) ; 
end 

jv= (Pprime-pl) . / (p2-pl) * ( jv6- jv5) + jv5; 
jl= (Pprime-pl) . / (p2-pl) * ( jl6- jl5) + jl5; 

else 

b 

a 

disp('hl and Pprime were not in the correct ranges for., 
f indjvjl .m' ) 
end 
end 
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